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Preface 



Einstein's general relativity remains, after almost a eentury, our best theory of gravitation. 
Its mathematical formalism and quantitative predictions will most likely play a crucial role 
also in future, especially on temporal and spatial scales that are relevant for astronomical 
purposes. This concerns mainly the description of stellar-mass compact objects as well 
as dark supermassive bodies in galactic nuclei. In spite of many evident successes, the 
general theory of relativity has not been tested yet with sufficient accuracy in the regime 
of strong gravity. Nevertheless, partial observational tests have been obtained, mainly 
by employing the methods of X-ray spectroscopy. Some questions have been answered, 
while others have emerged. In the near future, improved technology will allow to perform 
high-throughput high-resolution spectroscopy and polarimetry in X-rays with capabilities 
for detailed mapping of the close neighbourhood of compact objects. In this way, physical 
parameters of the central compact bodies will be measured and properties of accretion flows 
will be studied, including the innermost regions near the black-hole horizon. The aim of 
the present thesis is to contribute to the effort by developing an advanced computational 
tool for X-ray spectro-polarimetry in a strong gravity regime. 

We will turn our attention to nuclei of galaxies that harbour very massive objects, 
attracting a gaseous environment from their imminent surroundings. It appears beyond 
all reasonable doubt that these are black holes. The case for black holes residing in hearts 
of galaxies is not based on any single type of evidence. Instead, it is a mosaic of various 
observations that strongly supports the aforementioned interpretation. A similar situation 
arises on a vastly different scale of masses, namely in compact binary stars containing a 
stcUar-mass accreting black hole as one of its components. 

This thesis concentrates on a fragment of the whole picture. Our main goal is to 
discuss and improve the procedures and computational packages that are appropriate for 
spectroscopical studies of accreting black holes, including timing and polarimetry. The 
aim of this research is to discriminate between objects containing black holes from those 
containing other forms of compact bodies and to assess the possibility of determining 
physical parameters of these systems via the spectroscopy of accretion discs. We study the 
effects of curved space-time on light rays and specific spectral features which should arise 
when the gaseous environment emits radiation while falling onto the black-hole horizon. 
We thus develop an extended tool to fit X-ray data of real sources within a full general 
relativistic context. 
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Introduction 



It has been established observationally that significant amounts of glowing gas is often 
present in the cores of galaxies. This is particularly true in the case of very bright active 
galactic nuclei (AGN) and also in Galactic black-hole candidates (BHC). In spite of the 
big scale difference between the two types of cosmic objects, we can be rather confident 
that the accretion of gas proceeds in disc-like configurations in both cases. However, the 
sizes of the discs and their geometries are often quite uncertain. 



On black holes and accretion discs 



It is not so lon g ago that the general theory of relativity was brought to life by Albert 
EinsteinI ( 19ld ). Black holes belong to the first of its offspring, the static ones being 
described as early as in 1916 by Karl Schwarzschildl These objects are characterized by 
strong gravity and an imaginary horizon that enwraps a region from which not hing can 
escape , not even the light. Black holes have their ancestors in the ideas of John Michelll 
( 1784h and Pierre Simon Laplace (Ii79(tI : Il79flh . who defined similar objects more than one 
century before - stars with such large masses that the escape velocity from their surface is 
larger than the speed of light. But this would be a different story, a Newtonian one. Soon 
the family of black holes grew b igger, sever al new types were classified - black holes with 
a non-zero angular mo mentum (|KeriJ'l963l), with an ele ctric charge or wi th a hypothetic 



magnetic monopole LReissnei 1916. : .NordstronilEb lS: N ewman et al. 1965h 



However, electrically charged non-rotating black holes are not of direct astrophysical 
interest since, if the hole is immersed in a cloud of ionized or partially ionized plasma, 
which is a typical situation in cosmic environments and we assume it hereafter in this 
work, then any excess charge becomes rapidly neutralized by selective accretion. Never- 
theless, charged black holes have a role as a model for the real is tic bl a ck holes formed 
in a sp inning co l lapse, see Waldl (Il974 ): Blandford &: Znaiekl (Il977ll: Damour et al. 



(|l97Hh : IPamouT] (jl98(t : IPrice Thornel ()l9Rfih : iKaxas Vokrouhlickvl \\^^ . The 



propert ies of black holes have been , studi e d by numerous authors i n var i ous co r itexts ; 



Thorne. Price fc MacDoualdl (jl 986h : iKormendv Richstoriel (jl 99,^ : IWaJdl (jl 998h 



see e.g. 

for a survey of the problems and the general discussion of astrophysical connotations of 
black holes. As far as the Kerr (rotating) black holes are concerned, the possibility to 
extract its rotational energy either in the ergosphere by Penrose process ( Penrose fc Flovdl 
197lh or by Bl andford- Znaiek process, w hen the black hole is embedded in the electro- 
magnetic field ( Blandford fc Znaiekl 1977 ) , are really intriguing. 

The problem of the origin of black holes is very exciting as well. Oppenheimer consid- 
ered black holes as the last stage in the evolutio n of massive stars whe n after having burnt 
all their f uel neither the pressure of electrons dChandrasekhar 1931 ) nor the pressure of 
neutrons ( Oppenheimer fc Volkofj 19391 Ehoades fc R,uffinilll974h can stop them from col- 
lapsing. There are theories that the black holes may be nearly as old as the universe itself. 
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These primordial black holes ( Hawking 1971 ) may have been the seeds for the formation 
of galaxies. 

From the astrophysical point of view the motion of matter and the propagation of 
light near black holes has a special significance because the matter can emit light that can 
be observed and thus we can deduce what is going on out there. The emission f rom the 
stellar-mass and supermassive black holes themselves is rather low (jHawkinelll975ah . Many 
authors investigated these i ssues and calculated either the tirne-like and null geodesies in 
the b lack-hole space-times ( Shard 19791 : Chandrasekhar 19831: Bicak. Semerak &: Hadraval 
h 993h or even the more complex hydrodynamical models for the accreting matter or the 
outfl owing wind. The spherical accretion onto the black hole was firstlv considered as early 
I944I bv Bondi and Hovle. The disc-like accretion dLvnden-Bel] 19691: Pringle fc Eee? 



as m 



I972I : iNovikov fc Thorn jll973l: IShakura fc SunvaevI [19731 : IShapiro. Lightman fc Eardlev 



197(il ) is potentially much more efficient in converting the gravitational energy into the 



thermal energy which can be radiated away. It is also a more realistic configuration when 
we consider that the accreting matter has a non-zero angular momentum. 

The family of accretion discs is quite large. Let us just mention the Shakura-Sunyaev 
or standard discs (SS), the advection dominated accretion flow (ADAF), the convection 
dominated accretion flows (CDAF), the advection dominated inflows-outflows solution 
(ADIOS), the Shapiro-Lightman-Eardley discs (SLE) and the slim discs. Their behaviour 
depends mainly on the density, pressure, viscosity and opacity of the matter that they 
consist of as well as on the equa tion of state and their accretion rat e. For a basic review 
on acc retion disc theories see e.g. Kato. Fukue fc M ineshige ( 1998j) or Frank. King fc Eainel 
( 20021) . 

The interaction between an accretion disc and the matter above and below the disc 
may also play an important role in the accretion process and therefore the configuration of 
the black hole, accretion disc and corona seems to be more appropriate (fPaczviisk i 1978tj; 
Haardt fc MaraschilUB . 



As can be seen from the multitude of the disc models, some of the processes involved 
in the accretion are quite well understood. Nevertheless, there are still many open issues, 
one of them being the origin of viscosity in accretion discs. Turbulence and magnetic 
fields may play an important role here. Another important phenomenon is the existence 
of jets emerging from the vicinity of the central black hole, which may be due to the 
interaction between the black hole and the accretion disc. To tackle these problems much 
more complex three dimensional numerical models of ac cretion discs in curved space-time 
with magnetic field present have to be investigate d (e.g. Hawlev. Gammie fc Balbii3 199,4 
Nishikawa er^l200ll : Ide Villiers fc HawlevlEx)^ . 



If the disc is thick and/or dense, its own gravity, in particular farther away from the 
central black hole, cannot be neglected. These self-gravi tating di s cs we r e investig a ted b y 



various authors in various a p proximatio i is, star ting with Ostrikeil (Il964h. Bardeen 
Fishbone &: Moncried ( 1976[ ) , IPaczvhskil ( 1978ah and Karas. Lanza fc: Vokrouhlickv 



im), 



199a)- 



A very useful ^nd pertinent review of this subject was recently written by Karas, Hure 
Semerak L2004i ]. 



Linking theory with observation 

Last ce ntury w as rich in discovering highly energetic and violent objects in astronomy. 
In 1918 Curtisl discovered a long jet of material coming out from elliptica l galaxy M87 . 
Later, galaxies with particularly bright and compact nuclei were observed ()Sevfer i] ll943h 
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and they were classified as a new type of galaxies, the so-called Seyfert galaxies. But 
truly amazing discoveries came when radio astronomy fully developed. When quasars 



were observed for the first time in 1960s, they were called radio stars (jMatthews et al. 



196fll : Matthews fc Sandagelll963 ) because their emission in the radio part of the spectrum 



was en ormous and they appeared like point sources. It was three years later that Schmidt! 
(1963) realized that these objects were not hosted by our Galaxy but that they were rather 
distant extragalactic objects with highly redshifted spectral lines. Since then many strange 
radio sources with jets and lobes of various shapes and sizes have been discovered, mostly 
associated with optical extragalactic objects ~ centres of distant galaxies. The luminosity, 
size and variability of these sources indicate that highly energetic processes were once at 
work there. Because of their huge activity they are called active galactic nuclei. In 1960s 
it was realized that the only reasonable way of explaining the production of such large 
amounts of energy was to be sought in the gravitational energy conversi on by means of an 
accreting massive black hole ( Salpeteil ll964l : Zeldovich fc: Novikov 1964[ ). 

A new branch of astronomy has emerged, in which many novel discoveries were made: 
it was the X-ray astronomy and it developed quickly. The first observations were made by 
Geiger counters carried by rockets in 1950s and balloons in 1960s and later by specially 
designed X-ray astronomy satellites that have been orbiting the Earth si nce 1970s. The 



first Galactic X-ray source was discovered in 1962 - a neutron star Sco X-1 (iGiacconi et al. 



1962 ). Soon more discoveries followed. X-ray binary source Cygnus X-1 (|Bowver et al 



19651 ) being one of them. The fact that this object is a strong X-ray emitter and that the 



optical and X-ray emission varies on very short time scales (as short as one thousandth of a 
second) suggests that the companion might be a black hole. Many Galactic objects (most 
of them binaries containing neutron stars) and active galactic nuclei emitting X-rays have 
been observed since then. It is believed that AGN sources host a supermassive black hole 
with a mass of approximately 10^ — 10^ Mq and that some of the Galactic X-ray binaries 
contain a black hole with the mass equal to several solar masses. 

A breakthrough has come with the two most recent X-ray satellites, Chandra and 
XMM-Newton, which provide us with the most detailed images and spectra with an un- 
precedented resolution. Some of the observed objects exhibit redshifted, broad or narrow 
features in their spectra, thus suggesting that the emission is coming from the very close 
vicinity of a black hole. But still, the data acquired by the instruments on board these 
satellites are not as complete and detailed as would be necessary for an unambiguous in- 
terpretation and comparison with the theoretical models. New X-ray missions are being 
planned - Astro-E2 (already being assembled). Constellation- X and Xeus. Hopefully, 
with the data measured by the satellites involved in these observations, it will be possible 
not only to decide whether there is a black hole in the observed system but also to de- 
termine its properties, its mass and angular momentum, as well as the properties of the 
surrounding accretion disc. 



This thesis 

Nowadays both the theory describing accreting black-hole systems and the instruments for 
observing such objects are quite advanced. Computational tools for comparing acquired 
data with theories have been developed to such an advanced state that we are able now 
to find out properties of particular observed systems on the one hand and check if our 
models are good enough for describing these systems on the other. One of the tools 
used for processing X-ray data, which also implements the possibility to fit the data 
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within a certain set of models/ is xspec. This X-ray spectral fitting package contains 
various models for explaining the measured spectra but it lacks fully general relativistic 
models. One of the aims of this thesis has been to help to fill this gap and provide new 
general relativistic models that would be fast and flexible enough and would be able to fit 
parameters describing a black hole, mainly its mass and angular momentum. Although 
it is not possible to fit time resolved data within xspec and polarimetric data have not 
been available for X-ray sources yet (though they may become available when future X-ray 
missions materialize), we developed new models able to deal also with these issues. 

Our objective is to model X-ray spectra from accretion discs near compact objects. 
The observed spectrum depends not only on the local radiation emitted from the disc 
but it is also affected by strong gravitation on its way to the observer. Emission from 
different parts of the disc may be either amplified or reduced by gravitational effects. We 
will focus our attention on these effects in Chapter ^ where we deal with calculations of 
six functions that are needed for the integration of the total spectra of the accretion disc 
(including time varying spectra and polarimetric information). These functions include the 
Doppler and gravitational shifts, gravitational lensing, relative time delay, two emission 
angles and change of the polarization angle. 

The basic types of local emission coming from an accretion disc are summarized in 
Chapter 121 where also the equations for the observed fiux and polarization at infinity are 
derived. 

In Chapter 131 we describe new ky models that we have developed and that can be used 
either inside xspec or as standalone programs for studying the non-stationary emission 
from accretion discs or for studying polarimetry. These are models for a general relativistic 
line and Compton refiection continuum, general convolution models and a model for an 
orbiting spot. We also compare these models with some of the models already present in 

XSPEC. 

Some applications of the new models are summarized in Chapter |11 Firstly, we employ 
the models to high-quality X-ray data that are currently available for the Seyfert galaxy 
MCG-6-30-15. Then we calculate the flux from X-ray illuminated orbiting spot and polar- 
ization from an accretion disc illuminated from a primary source located above the black 
hole. 

There are several appendixes in this thesis. There is a summary of the basic equations 
for Kerr space-time in Appendix Appendix El describes the layout of data files used 
in this thesis. A detailed description of the integration routines that we have developed 
and that can be used inside the XSPEC framework for general relativistic computations is 
included in AppendixjO In Appendix^an atlas of contour figures of all transfer functions 
is shown. 

This thesis describes my original work except where references are given to the re- 
sults of other authors. Parts of the research were carried out in collaboration and pub- 
lished in papers, as indicated in the text. In pa rticular, various features and usa ge of 
the newlv deve l oped c ode have been described in Dovcia k. Karas fc Yaaoob* ('2004') and 
Dovciak et al. J2004bh. Chapte r IH i s based on papers IPo vciak. Karas &_Matt (2004); 



uovciaK et aLi iizuu4di). unapte r |4| i s based on papers LUo vciaK. ivaras &_Matt [zum) 
Dovciak. Karas fc YaaoobI ( 2004h and Dovciak et al\ (l2004al) . Another applic ation to ran- 



dom emitting spots in the accretion disc can be found in lCzernv et al. (|2004h but it is not 
included in the present thesis. 

^Notice that, in accordance with traditional terminology, we use the term 'model' in two slightly dif- 
ferent connotations: first, it can mean a physical model of an astronomical system or a process that 
is under discussion; other possible meaning of the word is more technical, referring to a computational 
representation, or a routine that is employed in order to compare actual observational data with the theory. 



Chapter 1 



Transfer functions for a thin disc 
in Kerr space-time 



Photons emitted from an accretion disc are affected by gravitation in several ways. They 
change their energy due to the gravitational and Doppler shifts because photons are emit- 
ted by matter moving close to the source of strong gravitation. The cross-section of light 
tube changes as the photons propagate in the curved space-time. This effect is strongest 
for high inclinations of the observer when trajectories of photons coming from behind the 
black hole are bent very much. Bending of light also influences the effective area from 
which photons arrive due to different emission angles. Here, aberration caused by motion 
of the disc matter plays its role as well. All three of these effects, the g'-factor, lensing and 
emission angle, influence the intensity of light that the observer at infinity measures. 

If the local emission from the disc is non-stationary one must take into account the 
relative delay with which photons from different parts of the disc arrive to the observer. 
A part of the local emission may be due to the reflection of the light incident on the disc 
from the corona above. In this case the disc radiation may be partially polarized. The 
polarization vector changes as the light propagates through the curved space-time. The 
change of polarization angle that occurs between the emission of light at the disc and its 
reception by the observer at infinity is needed for calculation of the overall polarization 
that the observer measures. When calculating the local polarization we need to know the 
geometry of incident and emitted light rays. Therefore we need to know also the azimuthal 
emission angle of the emitted photons. Thus six functions are necessary to transfer the 
local flux to the observer at infinity. The exact definitions for all of them will be given in 
next sections. For simplicity we will refer to them as transfer functions. 

Various authors have computed radiation from matter moving around a black hole in 
differe nt approximations and in different parameter space. De Felice, Nobili &: Calvani 
(jl974h computed the effects of gravitational dragging on the electr omagnetic radiation 
emitted by particles moving on bound orbits around a Kerr black hole. ICunning 



studied the combined effects of the gravitational and Doppler shifts together with the 
gravitational lensing effect on the X-ray radiation from an accretion disc in the strong 
gravitational field of a black hole. He introduced a concept of a transfer function where 
he includes all of the above-mentioned effects into a single function that describes the 
overall influence of the gravitational field on light rays emerging from the disc. He 
also studie d self-irradiation of the disc and its impact on the disc emission (Cunning- 
ham 1976| ). Since then numerous authors have investigated emission from disc (line, 
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blackbody, reflected, etc.) in strong gravity 


- Gerbal & Pelat (Il9811: Asaoka ( 


198C 


); 
); 
); 


Fabian et al\ (Il98c 


1: Koiima (Il991): Laod (Il9911: Karas. Vokrouhlickv k PolnarevI ( 


1992 


Viereutz (199.'^): 


Bao. Hadrava & 0stgaardl (1 


1994h: iHameurv. Marck k PelatI ( 


1994 


Zakharovl h9Mh: "Sneith. Riffert <^ Ruder] 


h99,4l: BroTTilev. Chen Miller' ( 


1997 


); 


Martocchia. Karas & Matt (.200(11: ICadez. Calvani k Fanton (|200.1): .Beckwith & Bond 



()2004l ). In the situation when full spectral resolution is not available, power spectra 
of light curves provide us with partial information on the emission from accretion discs 



(see e.g. Abramowicz et al. 199ll : Schnittman Bertschingei 2004). Other works also 



included the studies of n et po l arization from accretion discs , here we should m ention 
Connors. Piran fc Starlj (Il98nh: iLaor. Netzer fc PiranI (Il99(ih: IChen fc EardlevI ([l991'l : 



Matt] (1199.^): iMatt. V^hS^^ k Rosj (jl99.-^: lAgofl (|l997l): iBan et nl\ (jlMd); 

Qgura. Ohuo k Koiimal (|2000|). All of the above-mentioned authors used various methods 



for computing the transfer of photons from the accretion disc to the observer far away from 
the central black hole. Some of them assumed the black h ole to be non-rotating, others 
perfor med their calculations for a rotating K err black hole. iKaras. Lanza k Vokrouhlickv 
( 1995h and Usui. Nishida k Eriguchil (|l998f ) explored a substantially more complex case 
when self-gravity of the disc is not neglected and photons move in a space-time, the met- 
ric of which is itself derived from a numerical solution of Einstein's equations. Similar 
approach is appropriate also for ray-tracing in the field of a fast rotating neutron star, 
and in this way these techniques are pertinent to the study of sources of quasi-periodic 
oscillations in X-ray binaries. 

Our numerical codes for computing the emission from accretion discs (see Chapter ^ 
have to be fast, so that they are suitable for fitting data. That is the reason for which we 
have decided to pre-calculate the transfer functions and store them in the form of tables in 
a FITS^ file (see Appendix lB)) . We have chosen to compute the transfer of photons from the 
infinity (represented by setting the initial radius n = 10^^ in our numerical calculations) to 
the disc numerically by solving the equation of the geodesic. For this purpose the Bulirsch- 
Stoer method of integration has been use d (alternatively, the semi- analytical method of 
elliptical integrals could be used, see e.g. Ranch k BlandfordI 1994| ). An optically thick 
and geometrically thin flat disc has been assumed. Only null geodesies starting at the 
observer at infinity and ending at the equatorial plane of the black hole, without crossing 
this plane, has been taken into account (higher order images of the disc have not been 
computed). The Kerr space-time has been assumed and the medium between the disc 
and the observer has been supposed to be optically thin for the wavelengths that we are 
interested in (X-rays in our case). Special Kerr ingoing coordinates, which are non-singular 
on the horizon and bring spatial infinity to a finite value (to zero), have been used (see 
AppendixEJ)- This has enabled us to calculate the transfer functions with a high precision 
even very close to the horizon of the black hole. By using these coordinates the Boyer- 
Lindquist azimuthal coordinate (p has been unfolded to the Kerr coordinate (/9k , which 
is more precise when interpolating between the computed values on the disc. Some of 
the transfer functions depend on the motion of the matter in the disc, which has been 
assumed to be Keplerian above the marginally stable orbit and freely falling below 
it, where it has the same energy and angular momentum as the matter which is on the 
marginally stable orbit. All of the transfer functions except lensing have been actually 
calculated analytically (see formulae in the next sections) but with the numerical mapping 
of the impact parameters a and /3 to the disc coordinates r and (/Jk- The lensing has 



^Flexible Image Transport System. See e.g. lHanisch et ail Jij 

specifications. 



or |http ://fits.gsfc. nasa . go v/| for 
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CHAPTER 1. TRANSFER FUNCTIONS FOR A THIN DISC 



been computed by numerical integration of the equation of the geodesic deviation by the 
same method as described above. The initial conditions for the numerical integration are 
summarized in Appendix IA.4I In numerical computations, approximately 10^ geodesies 
have been integrated. They have covered the disc from the horizon to an outer radius of 
''out ~ 1000 in such a way that the region closest to the horizon has been most densely 
covered. Thus we have obtained values of transfer functions on a non-regular grid over 
the disc which have been interpolated to a regular grid in r and coordinates using 
Delaunay triangulation. Graphical representation of the transfer functions can be seen in 
Appendix IdI 

In the following sections we summarize the equations used for the evaluation of par- 
ticular transfer functions over the disc. In this chapter as well as everywhere else in the 
thesis we use units where GM, = c = 1 with M, being the mass of the central black hole. 



1.1 Gravitational and Doppler shift 



The four-momentum of the photons e mitted from the disc is given by (see e.g. 
and lMisner. Thorne fc Wheeleil[l973[ ) 



Pi 
Pt 



[a{l -a) + (r2 + a^){r'^ + - al)//\]/r'^ , 
Rsgn{(r2 + a'- alf - A[(/ - af + q^W'/^r^ 

- a + a(r2 + - a/)/A]/r2 , 



(1.1) 
(1.2) 
(1.3) 
(1.4) 



where a is the dimensionless angular momentum of the black hole (0 < a < 1). In 
eqs. p.l() - (|1.4|) we employed usual notation (see Appendix^ for the exact meaning of the 
quantities). 

The combined gravitational and Doppler shift ((/-factor) is defined as the ratio of the 
energy of a photon received by an observer at infinity to the local energy of the same 
photon when emitted from an accretion disc 



Pot 



1 



(1.5) 



Here and fe denote the frequency of the observed and emitted photons respectively, 
and is a four-velocity of the matter in the disc. 

Contour graphs of the (7- factor in Appendix ^ show that the effect of rotation of the 
black hole is visible only in its vicinity (r < 10). Near the horizon the 5- factor decreases 
down to zero due to the gravitational redshift, whereas far from the black hole the Doppler 
shift prevails. Very far from the black hole, where matter of the disc rotates slowly, the 
g-factor goes to unity. 

Preliminary considerations on the line emission (its energy shift and its width) can be 
based on the extremal values of the energy shift, 5+ and which photons experience 
when arriving at the observer's location from different parts of the disc. This is particularly 
relevant for some narr ow lines whose redshift can be determined more accurate lv than if 
the line is broad (e.g. Turner et al. 20021 : Guainazzl 20031 : Yaaoob et al. 2003f ). Careful 
discussion of g± can, in principle, circumvent the uncertainties which are introduced by 
the uncertain form of the intrinsic emissivity of the disc and yet still constrain some of 
parameters. Advantages of this t echnique were pointed out already bv . Cunningham (1975) 
and it was further developed by Pariev. Bromlev &: Miller ( 200l[ ). 



1.1. Gravitational and Doppler shift 
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Figure 1.1: Maximum and minimum values of energy shift (left) and emission angle (right) 
of photons originating from different radii in the disc. The dependence on the black-hole 
angular momentum, a, is shown for three different inclination angles, 0^ = 0.1°, 45° and 
85°. Each curve corresponds to a fixed value of a in the range (0,1), encoded by line 
colours. 



Fig. 11.11 shows the extremal values of the redshift factor for the observer inclinations 
9o = 0.1°, 45° and 85°. These were computed along circles with radius r in the disc plane. 
Corresponding contour lines of constant values of g± in the plane o versus r can be seen in 
Fig. 11.21 In other words, radiation is supposed to originate from radius r in the disc, but 
it experiences a different redshift depending on the polar angle. Contours of the redshift 
factor provide a very useful and straightforward technique to determine the position of 
a flare or a spot, provided that a narrow spectral line is produced and measured with 
sufficient accuracy. The most direct use of extremal (^-values would be if one were able 
to measure variations of the line profile from its lowest-energy excursion (for g_) to its 
highest-energy excursion (for g^), i.e. over a complete cycle. Even partial information can 
help to constrain models (for example, the count rate is expected to dominate the observed 
line at the time that g = g^, and so the high-energy peak is easier to detect). While it is 
very difficult to achieve sufficient precision on the highly shifted and damped red wing of 
a broad line, prospects for using narrow lines are indeed interesting, provided that they 
originate close enough to the black hole and that sufficient resolution is achieved both in 
the energy and time domains. 

Fig. 11.21 also makes it clear that it is possible in principle (but intricate in practice) to 
deduce the a-parameter value from spectra. It can be seen from the redshift factor that 
the dependence on a is rather small and it quickly becomes negligible if the light from the 
source is dominated by contributions from r > 10. 
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Figure 1.2: Contour lines of extremal energy shift factors in the plane of the rotation 
parameter a versus radius r for observer inchnation = 45° (left) and 85° (right). Max- 
imum 5f-factor is shown at the top, minimum ^-factor at the bottom. Radius is in units 
of GM/(P. The horizon is shown in black. The curve a(r)|j,=r^^ (dashed) is also plotted 
across the contour lines. Notice that, in the traditional disc-line scheme, no radiation is 
supposed to originate from radii below marginally stable orbit r^^g. If this is the case then 
one must assume that all photons originate outside the dashed curve and so the effect of 
frame-dragging is further reduced. 



1.2 Emission angle 

Here we examine the angle of emission with respect to the disc normal. We assume the 
local frame co-moving with the medium of the disc. The cosine of the local emission angle 
is 

He = COS de = — , (1.6) 

where n" = — e"^^ are components of the disc normal. 

One can see from contours of plotted in the disc plane (Appendix ^ and from 
graphs of its maximum and minimum values (Fig. II. 1|) that gravitation influences visibly 
only those light rays that pass close to the black hole on their way to the observer, and 
especially for large inclination angles. Due to the effect of light bending and aberration 
there exists a point on the disc where the emission angle is zero, i.e. the local direction of 
the emission is perpendicular to the disc surface. On the other hand, near the horizon only 
the light rays emitted almost parallel to the disc plane can reach the observer. Very far 
from the black hole, where matter rotates slowly, the emission angle is determined mainly 
by special-relativistic aberration and it gradually approaches the value of the observer 
inclination. 



1.3. Gravitational lensing 
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Figure 1.3: Same as in Fig. 11.11 but for lensing (left) and overall transfer function (right). 



1.3 Gravitational lensing 

We define lensing as the ratio of the cross-section dS{ of the light tube at infinity to the 
cross-section dS± of the same light tube at the disc: 

d^f 1 



(1.7) 



The four-vectors Y^i and 1^2 are transported along the geodesic according to the equation 
of the geodesic deviation from infinity where they are unit, space-like and perpendicular to 
each other and to the four-momentum of light (for the exact definition of these vectors very 
far from the black hole see ApDendix lA.4|) . In eq. (jl.Tj) we have denoted the magnitude of 
a four- vector by || || and scalar product of t wo four- vectors by < , >. The c ross section 
of light tube is constant for all observers (see Schneider. Ehlers Falcol 19921 ). 

Lensing can significantly amplify the emission from some parts of the disc (located 
behind the black hole from the point of view of the observer and in Kerr ingoing coor- 
dinates). This is true mainly for observers with large inclination angles (see figures in 
Appendix IdI and Fig. II. 

When one wants to estimate the total effect that gravitation has on the intensity of 
light coming from different parts of the disc, one has to take into account all the three 
effects - the 5- factor, the emission angle and the lensing. These effects combine into a 
single function defined by (see eq. 1)2. Ill) ) 



F 



2 1 



(1.^ 



We call it the overall transfer function. See contour graphs in Appendix ^ and plots in 
Fig. 11.31 for the total gravitational amplification of photon flux emitted from a Keplerian 
accretion disc. 
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Figure 1.4: Same as in Fig. 11.11 but for maximum and minimum relative time delay. 



1.4 Relative time delay 

The relative time delay At is the Boyer-Lindquist time which elapses between the emission 
of a photon from the disc and its reception by an observer (plus a certain constant so that 
the delay is finite close to the black hole but not too close). We integrate the equation 
of the geodesic in Kerr ingoing coordinates and thus we calculate the delay in the Kerr 
ingoing time coordinate Atx- The Boyer-Lindquist time can be obtained from Atx using 
the following equation: 



dt = dtK 



1 + 



2r 

(r — r+)(r — r_) 



dr , 



(1.9) 



with r-t = 1 lb y/l — a"^ being inner (— ) and outer (+) horizon of the black hole. Then we 
integrate the above equation and define the time delay as 



At 
At 



AtK 
AtK 



r + 



2 



In — + In [(r — r_|_)(r — r_)] 



r- 1 



+ 21n(r - 1 



for a < 1, (1.10) 
for a = 1. (1.11) 



There is a minus sign in front of the brackets because the direction of integration is from 
infinity to the disc. 

For contour graphs of the relative time delay see figures in Appendix El and for plots 
of maximum and minimum delay at a constant radius see Fig. 11.41 



1.5. Change of polarization angle 
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(i) Let three- vectors pa, ^a, n'j^ and /a be the mo- 
mentum of a photon, normal to the disc, projec- 
tion of the normal to the plane perpendicular to 
the momentum and a vector which is parallelly 
transported along the geodesic (as four- vector) , 




respectively; 

(ii) let ^ A be an angle between and /a; 

(iii) let the quantities in (i) and (ii) be evaluated at 
the disc for A = \ with respect to the local rest 
frame co-moving with the disc, and at infinity for 
A = 2 with respect to the stationary observer at 
the same light geodesic; 



(iv) then the change of polarization angle is defined 
as ^ = ^2- ^1- 



Figure 1.5: Definition of the change of polarization angle ^. 



1.5 Change of polarization angle 

Various physical effects can influence polarization of light as it propagates towards an 
observer. Here we examine only the influence of the gravitational field represented by the 
vacuum Kerr space-time. The change of the polarization angle is defined as the angle 
by which a vector parallelly transported along the light geodesic rotates with respect to 
some chosen frame. We define it in this way because in vacuum the polarization vector is 
parallelly transported along the light geodesic. This angle depends on the choice of the 

local frame at the disc and at infinity. See Fig. 11.51 for an exact definition. 

The change of polarization angle is (see Connors Sz Starj|l977 : Connors. Piran Sz StarkI 



198nh 

tan* = ^, (1.12) 

where 

X = —{a — asin6o)Ki — I3k2 , (1-13) 
Y = [a — asin6o)K2 — (3k,i , (1-14) 

with Ki and K2 being components of the complex constant of motion (see Walker &; 
Penrose FflTnl) 

Ki = arplf -rlapi-ir"^ + a^)pt]f -r{r'^ + a^)plf'P , (1.15) 
^^2 = -rplf + r[pi-apt]r + arpir. (1.16) 

Here the polarization vector is a four-vector corresponding to the three-vector /i from 
Fig. II. 51 which is chosen in such a way that it is a unit vector parallel with (i.e. ^'i = 0) 

For contour graphs of the change of the polarization angle see figures in Appendix IdI 
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1.6 Azimuthal emission angle 

We define the azimuthal emission angle as the angle between the projection of the three- 
momentum of the emitted photon into the equatorial plane (in the local rest frame co- 
moving with the disc) and the radial tetrad vector: 



spectral profiles, but it appears in discussions of polarimetry. This is also the reason why 
we have computed $e- Like the quantities discussed previously, the resulting values of the 
azimuthal angle depend on the adopted geometry of the source and the rotation law of 
the medium. They are determined by mutual interplay of special- and general-relativistic 
effects. Therefore, we remind the reader that our analysis applies to geometrically thin 
Keplerian discs residing in the equatorial plane, although generalization to more compli- 
cated situations should be fairly straightforward. For contour graphs of the change of the 
azimuthal emission angle see figures in Appendix IdI 





The explicit formula for this angle is not necessary for computations of light curves and 



Chapter 2 



Radiation of accretion discs in 
strong gravity 



There is now plausible evidence that the emission in some active galactic nuclei and some 
Galactic X-ray binary black-hole candidates originates, at least in part, from an accretion 
disc in a strong gravitational field. A lively debate is aimed at addressing the question 
of what the spectral line profiles and the associated continuum can tell us about the 
central black hole, and whether they can be used to constrain parameters of the accretion 
disc in a nearby zo ne, about ten gravit at ional radii or less from th e centre. For a recent 
review of AGNs see iFahian p.t fl/T|2n nn): Re vnolHs Nowa.kl (l2nn,'l) . and references cited 
therein. For BHCs see lMiller et al\ (|2002l ^ : iMcChntock RemillardI (|2003h and references 
therein. In several sources there is indication of iron Ka line emission from within the 
last stable orbit of a Schwarzschild black hole, e.g. in the Seyfert galaxy MCG-6-30-15 
fseejlwasawa, et rfi.lll994 iFahia.n et nJlbood ^^^^^^^ iMartorrhia. Matt KaTa^ 
20021 1 or from the region near abo ve the marginally stable orbit, a s in the case of the X-ray 



transient source XTE J1650-500 (jMimutti. Fabian fc Mi'lieill2004h that has been identified 
with a Galactic black-hole candidate. In other cases the emissi on appears to arise fart her 
from the b lack hole fe.g. in the microquasa r GRS 1915+105, see lMartocchia et fli.ll2n02l : for 
AGNs see Yaqoob fc Padmanabhaiil 20041 and references therein) . Often, the results from 
X-ray line spectroscopy are inconclusive, especially in the case of low spectral resolution 
data. For example, a spinni ng black hole is all owed but not required by the line model of 
the microquasar V4641 Sgr (|Miller et fli.ll2002h . The debate still remains open, but there 
are good prospects for future X-ray astronomy missions to be able to use the iron Ka line 
to probe the space-time in the vicinity of a black hole, and in particular to measure the 
angular momentum, or the spin, associated with the metric. 

One may also be able to study the 'plunge region' (about which very little is known), 
between the event horizon and the last stable orbit, and to deter mine if anv appreciable 
contr i bution to the i ron Ka line emission originates from there ( Revnolds Sz BegelmanI 



19971 : iKrolik Hawlev .2002 ). In addition to the Fe K lines, there is some evidence for 



relativis tic soft X-ray emis sion lines due to the Lya transitions of oxygen, nitrogen, and 
carbon ( Mason et al. 20o3). altho ugh the observational support for this interpretation is 
still controversial ()Lee et ai?.ll200lh . 

With the greatly enhanced spectral resolution and throughput of future X-ray astron- 
omy missions, the need arises for realistic theoretical models of the disc emission and 
computational tools that are powerful enough to deal with complex models and to allow 
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Model 






Effects that are taken into account 








Energy shift / 


a 


Non-axisymmetric 


Emission from 


Timing 


Polarization 




Lensing effect 




emission region 


plunge region 






DISKLINE 


yes/no 





no 


no 


no 


no 


LAOR 


yes/yes 


0.998 


no 


no 


no 


no 


KERRSPEC 


yes/yes 


(0,1)+ 


yes 


no 


no 


no 


KY 


yes/yes 


(0,1) 


yes* 


yes 


yes 


yes 



t The value of the dimensionless a parameter is kept frozen. 

* A one-dimensional version is available for the case of an axisymmetric disc. In this axisym- 
metric mode, KY still allows a and other relevant parameters to be fitted (in which case the 
computational speed of ky is then comparable to laor). The results can be more accurate than 
those obtained with other routines because of the ability to tune the grid resolution. 



Table 2.1: Basic features o f the new mo d el in c ompar i son w i th other black-hole disc- line 
models. Fo r references se e iFabian et al. (|l989h . Laoi ( 199l[ ). Martocchia. Karas fc Mat j 
( 2nn(t and bovciak et al\ l\20M\h . respectively. 



actual fitting of theoretical models to observational data. It is worth noting that some of 
the current data have been used to address the issue of distinguishing between different 
space-time metrics around a black hole, however, the current models available for fitting 
X-ray data are subject to various restrictions. 

In this chapter and the next one we describe a generalized schem e and a code which 
can be used with the standard X-ray spectral fitting package xspec ( ArnaudI 1996[ ). We 
have in mind general relativity models for black- hole accretion discs. Apart from a bet- 
ter numerical resolution, the principal innovations compared to the currently available 
schemes (see Tab. 12. Ij) are that the new model allows one to (i) fit for the black-hole 
spin, (ii) study the emission from the plunge region, and (iii) specify a more general form 
of emissivity as a function of the polar coordinates in the disc plane (both for the line 
and for the continuum). Furthermore, it is also possible to (iv) study time variability 
of the observed signal and (v) compute Stokes parameters of a polarized signal. Items 
(i)-(iii) are immediately applicable to current data and modelling, while the last two men- 
tioned features are still mainly of theoretical interest at the present time. Time-resolved 
analysis and polarimetry of accretion discs are directed towards future applications when 
the necessary resolution and the ability to do polarimetry are available in X-rays. Thus 
our code has the advantage that it can be u sed with tim e-resolve d data for reverb eration 
studies of relativistic accr etion discs fiiStellaiiil990 : Revn olds et aL,1999 : Ruszkowski 200cj; 



Govder fc Lasenbv 



2f)n4h . Also polarimetric analysis can be performed, and this will 
be extremel y useful because it can add ve r y specific information on strong-gravitationa l 
field effects (jConnors. Piran Starklll 98(11 : iMatt. Fahian RnH4Umi : hfio et »i.lll9"97h . 

Theoretical spectra with temporal and polarimetric information can be analysed with the 
current version of our code and such analysis should provide tighter constraints on future 
models than is currently possible. 



2.1 Basic spectral components of X-ray sources 

There are several components in the spectra of X-ray sources. Not all of them are always 
present and some of them are more prominent in certain objects than the others. 



2. 1 . Basic spectral components of X-ray sources 
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black ' hole accretion disc black i hole accretion disc 



Figure 2.1: Reflection models: a) lamp-post model; b) diffuse corona model. 



One of the characteristic spectral features that is almost always present is a power-law 
component. It is assumed that this feature results from inverse Compton scattering of 
thermal photons in a hot corona above the accretion disc. Usually two different configura- 
tions are considered. The first one assumes a patch of hot corona placed on the rotational 
axis of the central black hole at some height above it (lamp-post model). It is usually sup- 
posed to be an isotropic and point-like source of stationary primary power-law emission. 
The second configuration is a diffuse optically thin hot corona near above the accretion 
disc. For simple sketches of both configurations see Fig. 12.11 In both cases the shape 
of this primary power-law continuum (in the rest frame of the corona) is not affected by 
the relativistic effects acting on photons during their journey to the observer at infinity. 
These effects change only the normalization of the spectra. This component extends up 
to a cut-off energy of some tenths or a few hundreds keV. 

Often, mainly in the spectra of active galactic nuclei, additional continuum emission, 
a "hump", is added to the primary component. It is assumed that this is due to the 
reflection of the primary emission from the illuminated disc. The shape of the reflected 
continuum in the local rest frame co-moving with the accretion disc depends mainly on 
photoelectric absorption and Compton scattering of photons hitting the disc. The local 
emission is then smeared by relativistic effects - this concerns mainly its sharp features 
(e.g. iron edge). The shape of the observed spectra also depends on the illumination 
of the disc. This differs for the lamp-post model and the diffuse corona model. In the 
former case the radial dependence of illumination is determined by the height at which the 
patch of corona is placed. In the latter case the illumination of the disc depends on the 
emissivity of the diffuse corona near above the disc which may have quite a complicated 
radial dependence (but usually is assumed to be decreasing with radius as a power law). 

Spectral lines are an important feature observed in X-ray spectra. We assume that 
the origin of lines is the same as in the previous case ~ the illumination of the cold disc 
by primary (power-law) emission and reflection, in this case, by fluorescence. Originally 
narrow spectral lines are blurred by relativistic effects and thus they become broad, their 
width being as large as several keV in several sources. The most prominent examples are 
the iron lines Ka and K/?. 

In the present thesis we are concerned with the spectral components described above 
but here, we should not omit to mention other two important components - the thermal 
emission and the warm absorber. The thermal emission is more noticeable in the spectra of 
X-ray binaries with the black- hole candidate as the central object. This is due to the fact 
that in these sources the temperature of the accretion disc is much higher and therefore the 
black body emission extends as far as the soft X-ray energy band. This type of emission 
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a) b) c) 




Figure 2.2: Denomination of various elements of solid angles and areas defined in the text: 
a) the light source appears to the observer to be point-like; b) the hght rays received by 
the detector are coming from different parts of the disc (closer view of the disc than in 
previous figure) ; c) area of a light tube changes as the light rays travel close to the black 
hole (the disc is edge on). 



manifests itself in the spectra of active galactic nuclei as a "big blue bump" in the optical 
and UV frequencies. 

If there is ionized matter, a warm absorber, in the line of sight of the observer then 
the absorption spectral features are present in the observed spectra. These can include 
edges or resonant absorption lines due to oxygen and other ions. 

2.2 Photon flux from an accretion disc 

Properties of radiation are described in terms of photon numbers. The source appears as a 
point-like object for a distant observer, so that the observer measures the flux entering the 
solid angle df^o) which is associated with the detector area dSo^D^ dQo (see Fig. 12.2b ). 
This relation defines the distance D between the observer and the source. We denote the 
total photon flux received by a detector, 

^o{E)^^^= jd^N,{E/g)g\ (2.1) 

where 

(2.2) 

dr dS\ a\l\ 

is a local photon flux emitted from the surface of the disc, (in{E) is the number of photons 
with energy in the interval {E,E + dE) and g = E/E\ is the redshift factor. The local 
flux, Ni{Ei), may vary over the disc as well as in time, and it can also depend on the local 
emission angle. This dependency is emphasized explicitly only in the final formula (|2.11|) . 
otherwise it is omitted for brevity. 

The emission arriving at the detector within the solid angle dfi (see Fig. 12.2b ) originates 
from the proper area dSi on the disc (as measured in the rest frame co-moving with the 
disc). Hence, in our computations we want to integrate the flux contributions over a fine 
mesh on the disc surface. To achieve this aim, we adjust eq. ()2.1() to the form 

Here dS{ stands for an element of area perpendicular to light rays corresponding to the 
solid angle dO at a distance D, dS± is the proper area measured in the local frame of 



2.2. Photon Eux from an accretion disc 
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the disc and perpendicular to the rays, and d5 is the coordinate area for integration. We 
integrate in a two-dimensional slice of a four-dimensional space-time, which is specified 
by coordinates = Tr/2 and t = to — At with At being a time delay with which photons 
from different parts of the disc (that lies in the equatorial plane) arrive to the observer (at 
the same coordinate time to)- Therefore, let us define the coordinate area by (we employ 
coordinates t', 0, r, if with t' = t — At and At = At(r, 9,ip)) 



dS = \(fS/\ 



dt' ''^ 



|dV 



\9 



tfj,\ 



(2.4) 



We define the tensor d'^Sap by two four- vector elements dx^* = (dti,dr, 0, 0) and dx2 = 
(dt2, 0, 0, dy?) and by Levi-Civita tensor Sap-yS- The time components of these vectors, 
dti and dt2, are such that the vectors dx^ and dx2 lie in the tangent space to the above 
defined space-time slice. Then we obtain 



= W'steap dxf dx^l| = g''^-\\g^A^rdip = drdif.. 



(2.5) 



where g^y is the metric tensor and \\gni,\\ is the determinant of the metric. The proper 
area, dS±, perpendicular to the light ray can be expressed covariantly in the following 
way: 

i7-/dX/3 



dS^ 



(2.6) 



Here, dS± is the projection of an element of area, defined by d'^Sa/s, on a spatial slice 
of an observer with velocity C/" and perpendicular to light rays. C/" is four-velocity 
of an observer measuring the area dS±, and is four- momentum of the photon. The 
pro per area d5i corresponding to the same flux tube is identical for all observers 
(see Schneider. Ehlers Sz Falcolll992[ ). This means that the last equation holds true for 



any four- velocity C/", and we can express it as 



d^Sai3 + Pa dSs_ = , 



For a = t (note that d'^Str = d'^St^p 

dS± _ 
'dS' ~ 



- 0) we get 
1 dS^ 



g^^ d^S, 



te 



a 



Pt 



t, r, e, (f. 
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(2.7) 



(2i 



In the last equation we used the formula for the cosine of local emission angle fie, see 
eq. (|1.6|) . and the fact that we have chosen such an affine parameter of the light geodesic 
that Pt = —1. From eqs. ()2.3() . (|2.4() and 1)2. 8|) we get for the observed flux per unit solid 
angle 



dn{E) 



iVn 



dt d^o 

where A'o is a normalization constant and 

I = 



dr 



d(p Ni{E / g) g I He r, 



dSf 

dS^ 



(2.9) 



(2.10) 



is the lensing factor in the limit D — > oo (keeping D^dQ constant, see Fig. 12.2b ). 

For the line emission, the normalization constant Nq is chosen in such a way that the 
total flux from the disc is unity. In the case of a continuum model, the flux is normalized 
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to unity at a certain value of the observed energy (typically at -E = 1 keV, as in other 
XSPEC models). 

Finally, the integrated flux per energy bin, AE, is 



AN^{E,AE,t)= dEN^{E,t) = 

Je 

front j-(t>+^4> r-{E+l\E)/g 

= No dr dip dEiNi{Ei,r,ip,fie,t- At)g'^lficr, (2.11) 

Jr;„ J(h Je/q 



where At is the relative time delay with which photons arrive to the observer from different 
parts of the disc. The transfer functions g, I, /Xe and At are read from the FITS file 
KBHtablesNN . f its described in Appendix IB. II This equation is numerically integrated 
for a given local flux Ni{Ei, r,ip, ^(.,t — At) in all hereby described new general relativistic 
non-axisymmetric models. 

Let us assume that the local emission is stationary and the dependence on the axial 
coordinate is only through the prescribed dependence on the local emission angle /(/^e) 
(limb darkening/brightening law) together with an arbitrary radial dependence R{r), i.e. 

iVi(^i, r, if, fie,t - At) = Ni{Ei) R{r) f{p,). (2.12) 

The observed flux N^{E) is in this case given by 

/oo 
dE,N,{E,)G{E,E,), (2.13) 
-oo 

where 

rrovii r2iT 

G{E,Ei) = No drR{r) d^p f{fi,) I fi,r 6{E - gEi). (2.14) 

J Tin Jo 

In this case, the integrated flux can be expressed in the following way: 

f-E+AE 

AN^{E,AE) = / dEN^{E) = 

JE 

i-E+AE _ /Tout r2Tr _ fco 

= dENo drR{r) dip f {fi,) N,{E / g) g I fi, r dE,5{E, - E / g) = 

JE Jrin JO J -oo 

/Tout roo i-(E+AE)/Ei 

= No drR{r) dE,Ni{E{) dgF{g), (2.15) 

Jrin J-oo JE/Ei 

where we substituted g = E/E\ and 

/■27r 

F{g)= dipf{is,)g^lii,r5{g-g). (2.16) 
Jo 

Eq. 1)2.15(1 is numerically integrated in all axially symmetric models. The function dF{g) = 
dgF[g) has been pre-calculated for several limb darkening/brightening laws f{iJ.e) and 
stored in separate files, KBHlineNN . f its (see Appendix IB. 2|) . 
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2.3 Stokes parameters in a strong gravity regime 

For polarization studies, Stokes parameters are used. Let us define specific Stokes param- 
eters in the following way: 

Qu Vi/ . 

V^-^, qu = ^, Uy = —, v^ = ^^ (2-17) 

where 1^, Qu, Uy and Vy are Stokes parameters for light with frequency z^, ii^ is the energy 
of a photon at this frequency. Further on, we drop the index u but we will always consider 
these quantities for light of a given frequency. We can calculate the integrated specific 
Stokes parameters (per energy bin), i.e. Aio, Ago, Auq and Avo- These are the quantities 
that the observer determines from the local specific Stokes parameters i\, qi, ui and vi on 
the disc in the following way: 

dEiii{Ei)Fr, (2.18) 
dEi [qi{Ei) cos 2^' - ui{Ei) sin 2^'] Fr , (2.19) 



Aio{E,AE) 


= ""J 


f dS 


Aq^{E,AE) 


= N,j 


f dS 


Auo{E,AE) 


= N„j 


f dS 


Avo{E,AE) 


= N„j 


f dS 



dEivi{Ei)Fr. (2.21) 

Here, F = F{r,ip) = I fic is a transfer function, ^ is the angle by which a vector 
parallelly transported along the light geodesic rotates. We refer to this angle also as a 
change of the polarization angle, because the polarization vector is parallelly transported 
along light geodesies. See Fig. 11.51 for an exact definition of the angle ^. The integration 
boundaries are the same as in eq. (|2.11|) . As can be seen from the definition, the first 



specific Stokes parameter is equal to the photon flux, therefore, eqs. (|2.11|) and (|2.18j) are 
identical. The local speciflc Stokes parameters may depend on r, ^p, fig and t — At, which 
we did not state in the eqs. (|2.18|l - (|2.21j) explicitly for simplicity. 

The specific Stokes parameters that the observer measures may vary in time in the 
case when the local parameters also depend on time. In eqs. (|2.18)) - (|2.21|) we used a law of 
tr ansformation of the S tokes parameters by the rotation of axes (eqs. (1.185) and (1.186) 



m 



Chandrasekhailll96(t 



An alternative way for expressing polarization of light is by using the degree of polar- 
ization Pq and two polarization angles Xo and .^oi defined by 



= Vq^o +< + vl/io , (2.22) 
tan2xo = Uo/Qo, (2.23) 
sin 2^0 = Vo/Vq^ + uI + v^. (2.24) 



2.4 Local emission in lamp-post models 

The local emission from a disc is proportional to the incident illumination from a power- 
law primary source placed on the axis at height h above the black hole. To calculate the 
incident illumination we need to integrate the geodesies from the source to the disc. 
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The four-momentum of the incident photons which were emitted by a primary source 
and w hich are striki ng t he disc at radius r is fsee eas. (IA.8l) - (|A.ll|) with / = and 9 = it/2 
or also Carter 19681 and Misner. Thorne Sz Wheeler 19731 ) 



p\ = l + 2/r + 4/A, (2.25) 

pr = R'^^^[{r' + a')'-A{a' + qlp'/r\ (2.26) 

pf = qL/r\ (2.27) 

pf = 2a/(rA), (2.28) 

where = sin^^L (/i^+a^)^/AL — is Carter's constant of motion with Al = h? — 2h+a^ , 
and with the angle of emission 0l being the local angle under which the photon is emitted 
from a primary source (it is measured in the rest frame of the source) . We define this angle 

by tan^L = —Pi}/Pi}^ where j>|y ^ = ^l"^ ^"^^ Pi} — Pl ^lIl with ^^'^ ^hji being the 
four-momentum of emitted photons and the local tetrad connected with a primary source, 
respectively. The angle is 0° when the photon is emitted downwards and 180° if the photon 
is emitted upwards. 

We denoted the sign of the radial component of the momentum by Rggn- We have 
chosen such an affine parameter for the light geodesic that the conserved energy of the 
light is —pit = —pht = 1- The conserved angular momentum of incident photons is zero 
(/l = 0). 

The gravitational and Doppler shift of the photons striking the disc which were emitted 
by a primary source is 



5L = - = ^^ = -^. (2.29) 



l^L ~ PI.JJI ~ Ul 

Here v\ and vi^ denote the frequency of the incident and emitted photons, respectively and 
C/l is a four- velocity of the primary source with the only non-zero component = 
Cosine of the local incident angle is 

Hi = I cos (5i I = , (2.30) 

Pl U^' 

where n° = —^"g-^ is normal to the disc with respect to the observer co- moving with the 
matter in the disc. 

We further define the azimuthal incident angle as the angle between the projection of 
the three-momentum of the incident photon into the disc (in the local rest frame co-moving 
with the disc) and the radial tetrad vector, 




-1 PiaGf 

$i = -R^g, arccos I . 1 + ^ , (2.31) 

' ' . ^2 Pifj.U'^ 

where Rgg^ is positive if the incident photon travels outwards {p^"^^ > 0) and negative if it 

travels inwards {p^^^ < 0) in the local rest frame of the disc. 

In lamp-post models the emission of the disc will be proportional to the incident 
radiation Nf{Ei) which comes from a primary source 

Nf{Ei) = NS{Ej^)^. (2.32) 



2.4. Local emission in lamp-post models 
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Here N^{Ei^) = Nqi^ -^l ^ is an isotropic and stationary power-law emission from a primary 
source which is emitted into a soUd angle d^L and which illuminates local area d^i on the 
disc. The energy of the photon striking the disc (measured in the local frame co-moving 
with the disc) will be redshifted 

^i = 5L^L. (2.33) 



The ratio d^li^/dSi is 

dr^L dr^L dS sin OhdOh dap dS 

dSi dS dSi dr dip dSi 

where (see eqs. (|23|) and (|27|) ) 



(2.34) 



dS = drdip = Id^^/I = -/^ ^dS±= ^ dSi_ . (2.35) 

Pi Pi 

Here we used the same space-time slice as in the discussion above eq. 1)2. 4|) and thus the 
element d'^Sap is defined as before, see eq. ()2.5() . Note that here the area d5_L is defined 
by the incident flux tube as opposed to dS± in eq. ()2.8() where it was deflned by the 
emitted flux tube. The coordinate area d5 corresponds to the proper area d5_L which is 
perpendicular to the incident light ray (in the local rest frame co- moving with the disc). 
The corresponding proper area (measured in the same local frame) lying in the equatorial 
plane will be 

d5i = |d%f I = |ef,) eW-^dVl = \9ee^' d\g\ = 

-1/2 PiiJ. -1/2 f^L J e fn QR\ 

= -9ee = ggg' — gh dS± . (2.36) 

Pi Pi 

Here we have used eq. ()2.4() for the tetrad components of the element d'^Sap, eqs. 1)2. 6|) and 

It follows from eqs. (|2.32|) - (|2.36|) that the incident radiation will be again a power law 
with the same photon index T as in primary emission 

Nf{Ei) = Noi E-^ , (2.37) 

with the normalization factor 



M - N .r-i /i 2/t sin^Ld^L , . 



The emission of the disc due to illumination will be proportional to this factor. 



Chapter 3 

New models for XSPEC 



We have developed several general relativistic models for line emission and Compton re- 
flection continuum. T he line mo dels are supposed to be more accurate and versatile 
than the LAOR model iLaoiJ 199lh . and substantially faster than the kerrspec model 
( Martocchia. Karas k MatdboOOl). Several model s of intrinsic emissivity were employed, 
including the lamp-post model iMatt et fli]ll992h . Among other features, these models 
allow various parameters to be fitted such as the black-hole angular momentum, observer 
inclination, accretion disc size and some of the parameters characterizing disc emissivity 
and primary illumination properties. They also allow a change in the grid resolution and, 
hence, to control accuracy and computational speed. Furthermore, we developed very 
general convolution models. All these models are based on pre-calculated tables described 
in Chapter n and thus the geodesies do not need to be calculated each time one integrates 
the disc emission. These tables are calculated for the vacuum Kerr space-time and for a 
Keplerian co-rotating disc plus matter that is freely falling below the marginally stable 
orbit. The falling matter has the energy and angular momentum of the matter at the 
marginally stable orbit. It is possible to use different pre-calculated tables if they are 
stored in a specific FITS file (see Add endix IB . II for its detailed description). 

There are two types of new models. The first type of model integrates the local disc 
emission in both of the polar coordinates on the disc and thus enables one to choose 
non-axisymmetric area of integration (emission from spots or partially obscured discs). 
One can also choose the resolution of integration and thus control the precision and speed 
of the computation. The second type of model is axisymmetric - the axially dependent 
part of the emission from rings is pre-calculated and stored in a FITS file (the function 
dF{g) = dgF{g) from eq. (|2.16)) is integrated for different radii with the angular grid 
having 20 000 points). These models have less parameters that can be fitted and thus are 
less flexible even though more suited to the standard analysis approach. On the other 
hand they are fast because the emission is integrated only in one dimension (in the radial 
coordinate of the disc). It may be worth emphasizing that the assumption about axial 
symmetry concerns only the form of intrinsic emissivity of the disc, which cannot depend 
on the polar angle in this case, not the shape of individual light rays, which is always 
complicated near a rotating black hole. 

There are several parameters and switches that are common for all new models (see 
Tab. EH): 



a/M - specific angular momentum of the Kerr black hole in units of GM/c (M is the mass 
of the central black hole) , 
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0. 
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0. 
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Table 3.1: Common parameters for all models. 



parameter 


unit 


default value 


minimum value 


maximum value 


phi 


deg 


0. 


-180. 


180. 


dphi 


deg 


360. 


0. 


360. 


nrad 




200. 


1. 


10000. 


division 




1. 


0. 


1. 


nphi 




180. 


1. 


20000. 


smooth 




1. 


0. 


1. 


Stokes 




0. 


0. 


6. 



Table 3.2: Common parameters for non-axisymmetric models. 



theta_o - inclination of the observer in degrees, 

rin-rh ~ inner radius of the disc relative to the black-hole horizon in units of GM/c^, 
ms - switch for the marginally stable orbit, 

rout-rh - outer radius of the disc relative to the black-hole horizon in units of GM/c^, 
zshif t - overall redshift of the object, 

ntable - number of the FITS file with pre-calculated tables to be used. 

The inner and outer radii are given relative to the black-hole horizon and, therefore, their 
minimum value is zero. This becomes handy when one fits the a/M parameter, because 
the horizon of the black hole as well as the marginally stable orbit change with a/M, and 
so the lower limit for inner and outer disc edges cannot be set to constant values. The ms 
switch determines whether we also want to integrate emission below the marginally stable 
orbit. If its value is set to zero and the inner radius of the disc is below this orbit then 
the emission below the marginally stable orbit is taken into account, otherwise it is not. 

The ntable switch determines which of the pre-calculated tables should be used for 
intrinsic emissivity. In particular, ntable = for KBHtablesOO . f its (KBHlineOO . f its), 
ntable = 1 for KBHtablesOl .fits (KBHlineOl .fits), etc., corresponding to non-axisym- 
metric (axisymmetric) models. 

The following set of parameters is relevant only for non-axisymmetric models (see 
Tab. EH: 

phi - position angle of the axial sector of the disc in degrees, 
dphi - inner angle of the axial sector of the disc in degrees, 
nrad - radial resolution of the grid. 
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Figure 3.1: Segment of a disc from which emission comes (view from above) 
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photon flux array photar contains^ ^ 



i = I/E, where / is the first Stokes parameter (intensity) 
q = Q/E, where Q is the second Stokes parameter 
u = U/E, where U is the third Stokes parameter 
V = V/E, where V is the fourth Stokes parameter 
degree of polarization, P = y^g^ + + v'^/i 
angle x[deg] of polarization, tan 2% = u/q 



angle C[deg] of polarization, sin 2^ 



^ the photar array contains the values described in the table and multiplied 
by the width of the corresponding energy bin 
* i? is the energy of the observed photons 



Table 3.3: Definition of the Stokes parameter. 



division - switch for spacing of radial grid (0 - equidistant, 1 - exponential), 
nphi - axial resolution of the grid, 

smooth - switch for performing simple smoothing (0 - no, 1 - yes). 
Stokes - switch for computing polarization (see Tab. 13. 3|) . 

The phi and dphi parameters determine the axial sector of the disc from which emission 
comes (see Fig. 13. ip . The nrad and nphi parameters determine the grid for numerical 
integration. If the division switch is zero, the radial grid is equidistant; if it is equal to 
unity then the radial grid is exponential (i.e. more points closer to the black hole). 

If the smooth switch is set to unity then a simple smoothing is applied to the final 
spectrum. Here N^iE^) = [N^iE^-i] + 2N^{E^) + N^iE^+i)]/4. 

If the Stokes switch is different from zero, then the model also calculates polarization. 
Its value determines which of the Stokes parameters should be computed by xspec, i.e. 
what will be stored in the output array for the photon flux photar; see Tab. 13.31 (If 
Stokes 7^ then a new ascii data file stokes.dat is created in the current directory, 
where values of energy E together with all Stokes parameters i, q, u, v, P, x[deg] and 
^[deg] are stored, each in one column.) 

A realistic model of polarization has been currently implemented only in the kylIcr 
model (see Section [3.2. II belowl. In other models, a simple assumption is made - the local 



3.1. Models for a relativistic spectral line 
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emission is assumed to be linearly polarized in the direction perpendicular to the disc (i.e. 
qi = iy = Ni and ui = vi = 0). In all models (including kylIcr) there is always no final 
circular polarization (i.e. u = ^ = 0), which follows from the fact that the fourth local 
Stokes parameter is zero in each model. 



3.1 Models for a relativistic spectral line 

Three general relativistic line models are included in the new set of xspec routines - non- 
axisymmetric Gaussian line model kygIline, axisymmetric Gaussian line model kygline 
and fluorescent lamp-post line model kyfIll. 

3.1.1 Non-axisymmetric Gaussian line model KYGlLINE 

The kygIline model computes the integrated flux from the disc according to eq. (|2.11|) . 
It assumes that the local emission from the disc is 



W) = -^/(/^e) exp 



\ 2 

E] — Erest \ 
1000 * 



\/2sigma / 
Ni{Ei) = jump ^^^^*^-^ip^^ _L_ /(^^) exp 



E] — Erest \ ^ 

1000 

V2 Sigma / 



for r > Tb , (3.1) 
for r < rb . (3.2) 



The local emission is assumed to be a Gaussian line with its peak flux depending on the 
radius as a broken power law. The line is defined by nine points equally spaced with the 
central point at its maximum. The normalization constant Nq in eq. 1)2.11(1 is such that 
the total integrated flux of the line is unity. The parameters defining the Gaussian line 
are (see Tab. 13.4(1 : 

Erest - rest energy of the line in keV, 
sigma - width of the line in eV, 

alpha - radial power-law index for the outer region, 
beta - radial power-law index for the inner region, 

rb - parameter defining the border between regions with different power-law indices, 
jump - ratio between flux in the inner and outer regions at the border radius, 
limb - switch for different limb darkening/brightening laws. 

There are two regions with different power-law dependences with indices alpha and beta. 
The power law changes at the border radius where the local emissivity does not need 
to be continuous (for jump 7^ 1). The rb parameter defines this radius in the following 
way: 

rb = rb X for rb > , (3-3) 

rb = -rb + rh for rb < , (3.4) 



where r^g is the radius of the marginally stable orbit and rh is the radius of the horizon 
of the black hole. 
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Table 3.4: Parameters of the non-axisymmetric Gaussian line model kygIline. Model 
parameters that are not common for all non-axisymmetric models are denoted by asterisk. 



The function /(/ie) = /(cos(5e) in eqs. (|3.1jl and l|3.2|) describes the limb darkening/ 
brightening law, i.e. the dependence of the local emission on the local emission angle. 
Several limb darkening/brightening laws are implemented: 

/(/Zc) = 1 for limb = 0, (3.5) 

/(^J = 1 + 2.06/Xe for limb = -1 , (3.6) 

/(/xc) = ln(l + /i-i) for limb = -2, (3.7) 

/(^e) = /ie""'' for limb/ 0,-1, -2. (3.8) 

Eq. (|3.5|) corresponds to the isotropic local emission, eq. (|3.6|) corresponds to limb 
darkening in an optically thick electron scattering atmosphere (used by Laor, see 
Phillips fc MeszarosI Il983 iLaor. Netzer fc Piranlll99nl : FLaon.lQQUl . and eq. ((HIZ|) corre- 



sponds to limb brightening predicted by some models of a fluorescent line emitted by a n 
accretion disc due to X-ray ihumination (jHaardl] l l993l : iGhisellini. Haardt~ Mat,tl[T99i. 

There is also a similar model kyg2line present among the new XSPEC models, which 
is useful when fitting two general relativistic lines simultaneously. The parameters are 
the same as in the kygIline model except that there are two sets of those parameters 
describing the local Gaussian line emission. There is one more parameter present, ratio21, 
which is the ratio of the maximum of the second local line to the maximum of the first 
local line. Polarization computations are not included in this model. 

Results of an elementary code test are shown in Figs. I3.2H3.4I The intrinsic emissivity 
was assumed to be a narrow Gaussian line (width (t=0.42 FWHM = 2 eV) with the am- 
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Figure 3.2: Comparative examples of simple line profiles, showing a theoretical line 
{Ej-est = 6.4 keV) with relativistic effects originating from a black-hole accretion disc. 
Different sizes of the annular region (axially symmetric) have been considered, assuming 
that the intrinsic emissivity obeys a power law in the radial direction (a = 3). Resolution 
of the line-emitting region was UyXn^ = 3000 x 1500 with a non-equidistant layout of 
the grid in Kerr ingoing coordinates, as described in the text. Left: Dependence on the 
inner edge. Values of Hn are indicated in the plot (the outer edge has been fixed at the 
maximum radius covered by our tables, rout = 10^). Right: Dependence on rout (with the 
inner edge at horizon, = r^). Other key parameters are: 9^ = 45°, a = 1.0. Locally 
isotropic emission was assumed in the disc co-rotating frame. 




Energy [keV] Energy [keV] 

Figure 3.3: More calculated line profiles, as in the previous figure. Left: Line profiles 
for different values of a. Notice the enhanced red tail of the line when the intrinsic 
emission is concentrated to the centre of the disc. Right: Line emission originating from 
four different azimuthal segments of the disc. This plot can serve as a toy model of non- 
axisymmetric emissivity or to illustrate the expected effects of disc obscuration. Obviously, 
the receding segment of the disc contributes mainly to the low-energy tail of the line 
while the approaching segment constitutes the prominent high-energy peak. These two 
plots illustrate a mutual interplay between the effect of changing a and the impact of 
obscuration, which complicates interpretation of time-averaged spectra. The radial range 
is rji < r < 10^ in both panels. 
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Figure 3.4: Dependence on the observer inclination, as given in the plot. Left: Non- 
rotating black hole, a = 0. Right: Maximally rotating black hole, a = 1. Other parameters 
as in Fig. 13.21 and 13.31 



plitude decreasing oc r~" in the local frame co-moving with the disc medium. Typically, 
the slope of the low-energy wing is rather sensitive to the radial dependence of emissivity. 
No background continuum is included here, so thes e lines can be compared with similar 
pure disc-line profiles obtained in previous papers (|Koiim ^ Il99ll: lLaoilll99lh which also 



imposed the assumption of axially symmetric and steady emission from an irradiated thin 
disc. Again, the intrinsic width of the line is assumed to be much less than the effects of 
broadening due to bulk Keplerian motion and the central gravitational field. 

Furthermore, Figs. I3.5H3.6I compare model spectra of widely used XSPEC models. In 
these examples one can see that the LAOR model gives zero contribution at energy be- 
low 0.1-Erest (in the disc local frame). This is because its grid has only 35 radial points 
distributed in the whole range 1.23<r<400. That is also why, in spite of a very efficient 
interpolation and smoothing of the final spectrum, the LAOR model does not accurately 
reproduce the line originating from a narrow ring. Also, the dependence on the limb dark- 
ening/brightening cannot be examined with this model, because the form of directionality 
of the intrinsic emission is hard-wired in the code, together with the position of the inner 
edge at r>r^s- This affects especially the spectrum originating near the black hole, where 
radiation is expected to be very anisotropic and flow lines non-circular. The diskline 
model has also been frequently used in the context of spectral fitting, assuming a disc 
around a non-rotating black hole. This model is analytical, and so it has clear advantages 
in XSPEC. Notice, however, that photons move along straight lines in this model and that 
the lensing effect is neglected. 

3.1.2 Axisymmetric Gaussian line model KYGLINE 

This model uses eq. 1)2. ISp for computing the disc emission with local flux being 

Ni{Ei) = 5{Ei-Erest), (3.9) 
R{r) = r-^^P""^. (3.10) 



The function dF{g) = dgF{g) in eq. (|2.16j) was pre-calculated for three different limb 
darkening/brightening laws (see eqs. H3.5|) ~ H3.7|) ) and stored in corresponding FITS files 
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Figure 3.5: Comparison of the output from xspec models for the disc-hne problem: LAOR, 
DISKLINE, and kygIline (line 1 corresponds to the same limb darkening law and a = 
0.9982 as in LAOR; line 2 assumes locally isotropic emission and a = as in diskline). 
Left panel: 6^0 = 30°, = 6, rout = 7. Right panel: 60 = 70°, = 100, rout = 200. 
Radial decay of intrinsic emissivity follows a = 3 power law. 




Figure 3.6: Similar to previous figure but in logarithmic scale and for three choices of 
the darkening law in kygIline - (1) /(/Uc) = 1 + 2.06/ie; (2) /(^e) = 1; (3) fifJ-e) = 
log(l + l/fJ-e)- Left panel: 60 = 30°; Right panel: Oq = 70°. In both panels, nn = rmsj 
rout = 400, a = 0.9982. 



KBHlineOO .fits - KBHline02 .fits. The local emission is a delta function with its maxi- 
mum depending on the radius as a power law with index alpha and also depending on the 
local emission angle. The normalization constant A^o in eq. (|2.15|) is such that the total 
integrated flux of the line is unity. 

There are less parameters defining the line in this model than in the previous one (see 
Tab.l!T31): 

Erest - rest energy of the line in keV, 
alpha - radial power-law index. 

Note that the limb darkening/brightening law can be chosen by means of the ntable 
switch. 

This model is much faster than the non-axisymmetric kygIline model. Although 
it is not possible to change the resolution grid on the disc, it is hardly needed because 
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Table 3.5: Parameters of the axisymmetric Gaussian line model kygline. Model param- 
eters that are not common for all axisymmetric models are denoted by asterisk. 



the resolution is set to be very large, corresponding to nrad = 500, division = 1 and 
nphi = 20 000 in the KYGlLiNE model, which is more than sufficient in most cases. (These 
values apply if the maximum range of radii is selected, i.e. rin=0, ms=0 and rout=999; 
in case of a smaller range the number of points decreases accordingly.) This means that 
the resolution of the kygline model is much higher than what can be achieved with the 
laor model, and the performance is still very good. 

3.1.3 Non-axisymmetric fluorescent lamp-post line model KYFiLL 

The line in this model is induced by the illumination of the disc from the primary power- 
law source located on the axis at height above the black hole. This model computes the 
final spectrum according to eq. ()2.11|) with the local photon flux 



^ ' ^ rdr y height2 + (a/M)2^ 



X exp 



E] — Erest\ ^ 

1000 

V 2 Sigma / 



(3.11) 



Here, gi^ is ratio of the energy of a photon received by the accretion disc to the energy 
of the same photon when emitted from a source on the axis, 0l is an angle under which 
the photon is emitted from the source (measured in the local frame of the source) and 
Hi = cos 6i is the cosine of the incident angle (measured in the local frame of the disc) - 
see Fig. 12.11 All of these functions depend on height above the black hole at which the 
source is located and on the rotational parameter a/M of the black hole. Values of gh, Oh 
and /ii for a given height and rotation are read from the lamp-post tables lamp. fits (see 
Appendix IB.3|) . At present, only tables for a/M = 0.9987492 (i.e. for the horizon of the 
black hole rh = 1.05) and height = 2,3, 4, 5, 6, 8, 10, 12, 15, 20, 30, 50, 75 and 100 are 
included in lamp. fits, therefore, the a/M parameter is used only for the negative values 
of height (see below). 

The factor in front of the function /(/Ui,/Lte) gives the radial dependence of the disc 
emissivity, which is different from the assumed broken power law in the kygIline model. 
For the derivation of this factor, which characterizes the illumination from a primary 
source on the axis see Section ITU 
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Table 3.6: Parameters of the fluorescent lamp-post line model kyfIll. Model parameters 
that are not common for all non-axisymmetric models are denoted by asterisk. 

The function /(//i,/Xc) is a coefficient of reflection. It depends on the incident and 
reflection angles. Although the normalization of this function also depends on the photon 
index of the power-law emission from a primary source, we do not need to take this into 
account because the final spectrum is always normalized to unity. Values of this function 
are read from a pre-calcul ated table which is stored in f luorescentJLine . f its file (see 
Matt. Perola fc Pirolll99ll and Appendix rOl). 

The local emission ()3.11|) is defined in nine points of local energy E\ that are equally 
spaced with the central point at its maximum. The normalization constant A^o in the 
formula 1)2. 11[) is such that the total integrated flux of the line is unity. The parameters 
defining local emission in this model are (see Tab. 13. 6|) : 

Pho Index - photon index of primary power-law illumination, 

height - height above the black hole where the primary source is located for height > 0, 
and radial power-law index for height < 0, 

Erest - rest energy of the line in keV, 

sigma - width of the line in eV. 

If positive, the height parameter works as a switch - the exact value present in the 
tables lamp. fits must be chosen. If the height parameter is negative, then this model 
assumes that the local emission is the same as in the KYGlLiNE model with the parameters 
alpha = —height, rb = and limb = —2 (Pholndex parameter is unused in this case). 

In Fig. 13. 71 we demonstrate that the broad iron emission lines due to illumination from 
the source placed on the axis depend heavily on the height where the "lamp" is located 
(left), as well as on the photon index of the primary emission (right). These graphs 
correspond to the iron Ka line with the rest energy of 6.4 keV. 



34 



CHAPTER 3. NEW MODELS FOR XSPEC 




0.1 1 10 0.1 1 10 

Energy [keV] Energy [keV] 

Figure 3.7: An example of a line profile originating from a disc in equatorial plane of a 
Kerr black hole (a = 0.9987, i.e. rh = 1.05) due to the illumination from a primary source 
on the axis. The kyfIll model was used. Left: Dependence of the line profile on the 
height (in GM/(?) of a primary source with photon index F = 2. Right: Dependence 
of the line profile on the photon index of the primary emission with a source at height 
3GM/c^ above the black hole. 



3.2 Compton reflection models 

We have developed two new relativistic continuum models - the lamp-post Compton 
reflection model kylIcr and the kyhIrefl model which is a relativistically blurred hrefl 
model that is already present in xspec. Both of these models are non-axisymmetric. 

3.2.1 Non-axisymmetric lamp-post Compton reflection model KYLlCR 

The emission in this model is induced by the illumination of the disc from the primary 
power-law source located on the axis at height above the black hole. As in every non- 
axisymmetric model, the observed spectrum is computed according to eq. ()2.11|) . The 
local emission is 



PhoIndex-lSin^Ld^L L 2 height 



A^l(Si)=5r ^ 1 1 - , . , °, ;,,,2 W,/^c) for height > 0, (3.12) 

rdr y height^ + (a/M) 

A?'i(^i) =r''"^sht J^Ei- fxe) for height < . (3.13) 

For the definition of g^, Ol and fii see Section lS. 1.31 and ADDendix lB.31 where pre-calculated 
tables of these functions in lamp, fits are described. 

The function f {E\; fic) gives the dependence of the locally emitted spectrum on 
the angle of incidence and the angle of emission, assuming a power-law illumination. 
This function depends on the photon index Pholndex of the power-law emission from a 
primary source. Values of this function for various photon indices of primary emission 
are read from the pre-calculated tables stored in ref spectra. fits (see Appendix IB. 5|) . 
These tables were calculat ed by the Monte Carlo simulations of Compton scattering in 
Matt. Perola fc Pirol (|l99l[ ). At present, tables for Pholndex = 1.5, 1.6, . . . , 2.9, 3.0 and 



for local energies in the range from 2keV to 300 keV are available. The normalization 
constant A'^o in eq. (|2.1H) is such that the final photon fiux at an energy of 3 keV is equal 
to unity, which is different from what is usual for continuum models in XSPEC (where the 
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Table 3.7: Parameters of the lamp-post Compton reflection model kylIcr. Model pa- 
rameters that are not common for all non-axisymmetric models are denoted by asterisk. 



photon flux is unity at 1 keV) . The choice adopted is due to the fact that current tables 
in ref spectra. fits do not extend below 2keV. 

The function /(£'i;//e), which is used for negative height, is an averaged function 
f{Ei;ni,fj,e) over fii 

fiEi; ^e) = [ d/xi fiEi; f,,, //e) • (3.14) 
Jo 

The local emission (|3.13|1 can be interpreted as emission induced by illumination from 
clouds localized near above the disc rather than from a primary source on the axis (see 
Fig. 12.11) . In this case photons strike the disc from all directions. 

For positive values of height the kylIcr model includes a physical model of polar- 
ization based on Rayleigh scattering in single scattering approximation. The specific local 
Stokes parameters describing local polarization of light are 





h + Ir 
{h + Ir) 




(3.15) 




h-Ir 
{h + Ir) 




(3.16) 


ni(^i) 


U 


Ni{Ei) , 


(3.17) 




= 0, 




(3.18) 



where the functions ly, Ir and U determine the angular dependence of the Stokes param- 
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Figure 3.8: General relativistic lamp-post Compton reflection model kylIcr with 
(dashed) and without (solid) iron lines Ka and K/3. The emission from the disc is induced 
by illumination from a primary source placed 2GM/(? (left) and 100 GM/c^ (right) above 
the black hole. 



eters in the following way 

h = + + 2(1 - - ^?) - ^^^e^ii^J{l - /^2)(i - m?) cos - 

-lil{l-lM^)cos[2{^,-^,)], (3.19) 

Ir = 1 + /uf + (1 - [2(«>e - $i)] , (3.20) 

U = -4/iiy(r^^l2^(r^^sin($e-^i)-2/ie(l-/U?)sin[2($e-$i)]. (3.21) 

Here and are the azimuthal emission and the incident angles in the local rest frame 
co-moving with the accretion disc (see Sections 11.61 and 12.41 for their defi nition). For the 
derivation of these formulae see the definitions (1.147) and eqs. (X.172) in IChandrasekhai 
(^196fli ). We have omitted a common multiplication factor, which would be cancelled any- 
way in eqs. ()3.15|) ~ H3.18() . The symbol ( ) in definitions of the local Stokes parameters 
means value averaged over the difference of the azimuthal angles <&c — We divide the 
parameters by {I\ + Ir) because the function /(Si; /i;, /ie), and thus also the local photon 
flux Ni{Ei), is averaged over the difference of the azimuthal angles. 

The parameters defining local emission in this model are (see Tab. I3.7]l : 

Pholndex - photon index of primary power-law illumination, 

height - height above the black hole where the primary source is located for height > 0, 
and radial power-law index for height < 0, 

line ~ switch whether to include the iron lines (0 - no, 1 - yes), 

E_cut - exponential cut-off energy of the primary source in keV. 

The tables ref spectra. fits for the function f{Ei;fj,i,^e) also contain the emission in 
the iron lines Ka and K/?. The two lines can be excluded from computations if the line 
switch is set to zero. The E_cut parameter sets the upper boundary in energies where the 
emission from a primary source ceases to follow a power-law dependence. If the E_cut 
parameter is lower than both the maximum energy of the considered dataset and the 
maximum energy in the tables for f{Ei; ^i,//e) in ref spectra, fits (300 keV), then this 
model is not valid. 
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Table 3.8: Parameters of the reflection kyhIrefl model. Model parameters that are not 
common for all non-axisymmetric models are denoted by asterisk. 



Examples of the Compton reflection emission component of the spectra with and with- 
out the fluorescent Ka and K/? lines are shown in Fig. 13.81 It can be seen that originally 
narrow lines can contribute substantially to the continuum component. 

3.2.2 Non-axisymmetric Compton reflection model KYHlREFL 

This model is based on an existing multiplicative hrefl model in combination with the 
POWERLAW model, both of which are present in xspec. Local emission in eq. ()2.11|) is 
the same as the spectrum given by the model hrefl*powerlaw with the parameters 
thetamin = and thetamax = 90 with a broken power-law radial dependence added: 

Ni{Ei) = r'^^P^^ HREFL*POWERLAW for r > ry, , (3.22) 

Ni{Ei) = jump r^^*^"^^P''^r"''^*^HREFL*POWERLAW for r < rt . (3.23) 

For a definition of the boundary radius rb by the rb parameter see eas. (I3.3p-(nni) . and 

aaoobi ( 



for a detailed description of the hrefl model see iDovciak. Karas Sz YaaoobI (|200J) and 
the XSPEC manual. The kyhIrefl model can be interpreted as a Compton-reflection 
model for which the source of primary irradiation is near above the disc, in contrast to 
the lamp-post scheme with the source on the axis (see Fig. 12. 1|) . The approximations 
for Compton reflection used in hrefl (and therefore also in kyhIrefl) are valid below 
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Figure 3.9: Comparison of the general relativistic kyhIrefl model with the non- 
relativistic hrefl(powerlaw). The relativistic blurring of the iron edge is clearly visible. 
The power-law index of the primary source is PhoIndex=2 (left) and PhoIndex=2.6 (right). 
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Figure 3.10: Comparison of the two new general relativistic Compton reflection models 
KYLlCR and kyhIrefl. The lamp-post kylIcr model is characterized by the height 
h above the disc where a primary source of emission is placed, the reflection kyhIrefl 
model is characterized by the radial power-law index a. Left: h = IGM/c^, a = 3.4. 
Right: h = lOOGM/c^, a = 1.5. 



~ 15 keV in the disc rest-frame. The normalization of the final spectrum in this model is 
the same as in other continuum models in xspec, i.e. photon flux is unity at the energy 
of 1 keV. 

The parameters defining the local emission in kyhIrefl (see Tab. 13.8]) are 

Pholndex - photon index of the primary power-law illumination, 
alpha - radial power-law index for the outer region, 
beta - radial power-law index for the inner region, 

rb - parameter defining the border between regions with different power-law indices, 
jump ~ ratio between fiux in the inner and outer regions at the border radius, 
Feabun - iron abundance relative to solar, 
FeKedge - iron K-edge energy, 

Escfrac - fraction of the direct flux from the power-law primary source seen by the 
observer. 



3.3. General relativistic convolution models 
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covf ac - normalization of the reflected continuum. 

Smearing of sharp features in continuum (iron edge) by relativistic effects is demon- 
strated in Fig. 13.91 where non-relativistic reflection model hrefl(powerlaw) is compared 
with our relativistic kyhIrefl model. Here, we set the radial power-law index a = 1 in 
kyhIrefl. Other parameters defining these models were set to their default values. 

We compare the two new relativistic reflection models kylIcr and kyhIrefl in 
Fig. 13.1(11 Note that the kylIcr model is valid only above approximately 2keV and the 
kyhIrefl model only below approximately 15keV. 

3.3 General relativistic convolution models 

We have also produced two convolution- type models, kyIconv and kyconv, which can be 
applied to any existing xspec model for the intrinsic X-ray emission from a disc around 
a Kerr black hole. We must stress that these models are substantially more powerful 
than the usual convolution models in xspec (these are commonly defined in terms of 
one-dimensional integration over energy bins). Despite the fact that our convolution mod- 
els still use the standard xspec syntax in evaluating the observed spectrum (e.g. ky- 
CONv(powERLAw)), our code accomplishes a more complex operation. It still performs 
ray-tracing across the disc surface so that the intrinsic model contributions are integrated 
from different radii and azimuths on the disc. 

There are several restrictions that arise from the fact that we use existing xspec 
models: 

- by local xspec models only the energy dependence of the photon flux can be defined, 

- only a certain type of radial dependence of the local photon flux can be imposed - 
we have chosen to use a broken power-law radial dependence, 

- there is no azimuthal dependence of the local photon flux, except through limb 
darkening law, 

- local flux depends on the binning of the data because it is defined in the centre of 
each bin, a large number of bins is needed for highly varying local flux. 

For emissivities that cannot be defined by existing XSPEC models, or where the lim- 
itations mentioned above are too restrictive, one has to add a new user-defined model 
to XSPEC (by adding a new subroutine to xspec). This method is more flexible and 
faster than convolution models (especially when compared with non-axisymmetric one), 
and hence it is recommended even for cases when these prefabricated models could be 
used. In any new model for XSPEC one can use the common ray-tracing driver for rel- 
ativistic smearing of the local emission: ide for non-axisymmetric models and idre for 
axisymmetric ones. For a detailed description see Appendixes \C1\ and l( ' . 21 

3.3.1 Non-axisymmetric convolution model KYdONV 

The local emission in this model is computed according to the eq. 1)2.11(1 with the local 
emissivity equal to 

Ni{Ey) = r-^^P''V(/"c) MODEL for r>rb, (3.24) 

iVi(^l) = jump r^"'^"^'P''V-^"*V(Aic) MODEL for r < . (3.25) 
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Table 3.9: Parameters of the non-axisymmetric convolution model kycIonv. Model 
parameters that are not common for all non-axisymmetric models are denoted by asterisk. 



For a definition of the boundary radius rb by the rb parameter see eqs. p.3() - (|3.4() and for 
the definition of different limb darkening laws f{ne) see eqs. (|3.5|) - (|3.8|) . The local emission 
is given by the model in the centre of energy bins used in xspec with the broken power- 
law radial dependence and limb darkening law added. Apart from the parameters of the 
MODEL, the local emission is defined also by the following parameters (see Tab. 13. 9|) : 

normal ~ switch for the normalization of the final spectrum, 
= - total flux is unity (usually used for the line), 

> - flux is unity at the energy = normal keV (usually used for the continuum), 

< - flux is not normalized, 
ne_loc - number of points in the energy grid where the local photon flux is defined, 
alpha - radial power-law index for the outer region, 
beta - radial power-law index for the inner region, 

rb - parameter defining the border between regions with different power-law indices, 
jump - ratio between the flux in the inner and outer regions at the border radius, 
limb - switch for different limb darkening/brightening laws. 

The local emission in each ky model has to by defined either on equidistant or exponential 
(i.e. equidistant in logarithmic scale) energy grid. Because the energy grid used in the 
convolution model depends on the binning of the data, which may be arbitrary, the flux 
has to be rebinned. It is always rebinned into an exponentially spaced energy grid in KY 
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Table 3.10: Parameters of the axisymmetric convolution model kyconv. Model parame- 
ters that are not common for all axisymmetric models are denoted by asterisk. 



convolution models. The ne_loc parameter defines the number of points in which the 
rebinned flux will be defined. 

3.3.2 Axisymmetric convolution model KYCONV 

The local emission in this model is computed according to eq. (|2.15|) with the local emis- 
sivity equal to 

Ni{Ei) = MODEL, (3.26) 
R(r) = r-^^P^^. (3.27) 

Except for the parameters of the model, the local emission is defined also by the following 
parameters (see Tab. 13.1(H) : 

alpha - radial power-law index, 

ne_loc - number of points in energy grid where local photon flux is defined, 

normal - switch for the normalization of the final spectrum, 
= - total flux is unity (usually used for the line), 

> - flux is unity at the energy = normal keV (usually used for the continuum), 
< - flux is not normalized. 

Note that the limb darkening/brightening law can be chosen through the ntable switch. 
This model is much faster than the non-axisymmetric convolution model KYClONV. 

3.4 Non- stationary model KYSPOT 

Emission in this model originates from a localized spot on the disc. The spot may either 
move along a stable circular orbit or fall from the vicinity of the marginally stable orbit 
down to the horizon. In the former case the spot is orbiting with a Keplerian velocity, 
in the latter it falls with energy and angular momentum of the matter on the marginally 
stable orbit. Because the emission changes in time this model cannot be included into 
XSPEC which cannot handle non-stationary problems. In spite of this we may include 
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Figure 3.11: Dynamical profile of an iron line (-Erest = 6.4keV) produced by an orbiting 
spot. Energy is on the ordinate, time on the abscissa. The horizontal range spans the 
interval of one orbital period. Time zero corresponds to the time when the observer 
receives photons from the spot at the closest approach to the observer. Left: Large spot 
(/? = 0.1) at the radius rgpot = SGM/c^. Right: Obscured small spot (/? = 4) at the 
radius rgpot = 5GM/c'^. Obscuration occurs between 315° — 540° measured from the 
closest approach. In both cases the dimensionless angular momentum parameter of the 
black hole is a = 1 and observer inclination is = 45°. The observed photon flux is 
colour-coded (logarithmic scale with arbitrary units). 



integrated emission received by the observer in a certain time span (corresponding to a 
certain part of the orbit) in future. 

In this model it is assumed that the corona above the disc is heated by a flare. Thus a 
hot cloud is formed, which illuminates the disc below it by X-rays. The disc reflects this 
radiation by Compton scattering and by fluorescence. In the current model we consider 
only the fluorescent part of the reflection. 

The local emission from the disc is 



iVi(Ei) = /(/ie) exp 



\ 2" 

El — Erest \ 



exp [-beta(Ar)2] for beta(Ar)^ < 4, (3.28) 



V2a , 

iVi(^l) = for beta(Ar)2>4. (3.29) 

Here the function /(/Uc) describes the limb darkening/brightening law (see Section 13.1.11 
for more details), Erest is the local energy of the fluorescent line, a = 2eV is its width, 
beta determines the size of the spot and (Ar)^ = r'^ + r^p^^ — 2 r r^pot cos {ip — ^^spot) with 
Tspot and (/?spot being polar coordinates of the centre of the spot. The emission is largest 
at the centre of the spot and decreases towards its edge as is obvious from eq. (|3.28|) . The 
Gaussian line in energies is defined by nine points equally spaced with the central point 
at its maximum. 

The local emission is defined by the following parameters (see Tab. I3.1T|) : 
Erest - rest energy of the line in keV, 

sw - switch for choosing the type of spot (1 - orbiting, 2 - falling), 
beta - parameter defining the size of the spot, 

rsp ~ radius (in GM/(?) at which the spot is orbiting or the initial radius in units of the 
marginally stable orbit if spot is in-falling. 



3.4. Non-stationary model KYSPOT 



43 



n T P TTl At AT* 


1 1 ml" 
Hill L 


LltldLllL Vdilut; 


iTi ml m 1 1 m t rpi 1 1 1 
llllllllllLllll Vdlut: 


m5i"vimnm Tmlnp 

llldJ^llllLllll VdlLlC 


a /M 

a./ n 


r' M In 


n QQ89 


n 
u. 


"1 


Llic Lcl_U 


deg 


ou. 


n 

u. 


Oa. 


3riii~2rli 


KjrlVl j C 


n 
u. 


n 
u. 


QQQ 


rout ~2rli 


n 1\/T 7^.2 


zu. 


n 
u. 


QQQ 

yyy. 


phi 


deg 


n 
u. 


-low. 


1 sn 


dptii 


deg 


OOU. 


n 
u. 


OOU. 


niTcLci 




OUU. 


"1 


1 nnnn 

iUUUU. 


QXVXS XUll 




n 
u. 


n 
u. 


"1 


nptii 




( ou. 


"1 

1. 


zuuuu. 


smooth. 




"1 


n 
u. 


"1 

1 . 


zshif "t 




n 
u. 


n QQQ 

-u.yyy 


1 n 

iU. 


n i" a V^n o 

11 L dU J. c 




n 


n 


QQ 


Ejx esx. 


KC V 




"1 


QQ 

yy. 


*sw 




"1 

1. 


"1 


9 

z. 


DeXEL 




n 1 

U. i 


n nni 

U.UUi 


iUUU. 


^2rsp 


n 1\/T 7^.2 


Q 
O. 


n 
u. 


1 nnn 

iUUU. 


*DSD 


de£ 


90. 


-360. 


360. 


*Norbits 




1. 


0. 


10 


*nt 




500. 


1. 


le6 


*limb 




0. 


-10. 


10. 


*polar 




0. 


0. 


1. 



Table 3.11: Parameters of the non-stationary model kyspot. Model parameters that are 
not common for all non-axisymmetric models are denoted by asterisk. 



psp - azimuthal angle in degrees where the spot is initially located (90° for the closest 
approach) , 

Norbits - number of orbits for an in-falling spot (not used for an orbiting spot), 
nt - time resolution of the grid, 

limb - switch for different limb darkening/brightening laws, 

polar - switch for polarization calculations (0 - not performed, 1 ~ performed). 

The KYSPOT model creates an ascii file kyspot.dat, where the dependence of the 
observed spectrum (columns) on time (rows) are stored. If the polar switch is set to 
unity then two other ascii files are created - kyspot_poldeg.dat and kyspot_psi.dat. 
Here, the dependence of the degree of polarization and the angle of polarization on time 
are stored. The polarization calculations in this model are based on the same assumptions 
as in the kygIline model. 

Examples of dynamical spectra of orbiting and in-falling spots can be seen in Figs. 13.111 - 
13.121 The first figure represents an evolving spectral line in a dynamical diagram for an 
orbiting spot of two different sizes one of which is being obscured on part of its orbit around 
the black hole. The characteristic radius of the larger spot is ~ 6.3GM/c^ and the radius 
of the smaller spot is ~ IGM/c^. The second fi gure represents the dynamical spectra 
of the spot falling from the marginally stable orbit down towards the black-hole horizon. 
The characteristic radius of the spot is ~ 0.2GM/c^. This figure reflects the effect of 
changing the energy shift along the spot trajectory. The in-spiral motion starts just below 
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Figure 3.12: Dynamical profile of an iron line produced by an in-spiralling small spot 
(/3 = 100) in free fall. Left: Trajectory (in the equatorial plane) of the spot, defined by a 
constant energy and angular momentum during the in-fall. The observer is located to the 
top of the figure. Right: Spectrum of the spot. Energy is on the ordinate, time on the 
abscissa. The horizontal range spans the interval of 2.4 orbital periods at the corresponding 
initial radius, rgpot = 0.93 (a). Here, the dimensionless angular momentum parameter 
of the black hole is a = 0.9, observer inclination is 6^ = 45°. The observed photon flux is 
colour-coded (logarithmic scale with arbitrary units). Gradual decay of the signal and an 
increasing centroid redshift can be observed as the spot completes over two full revolutions 
and eventually plunges into the black hole. 



the marginally stable orbit (rsp = 0.93) and it proceeds down to the horizon, maintaining 
the specific energy and angular momentum of its initial circular orbit at r^s(a). One can 
easily recognize that variations in the redshift could hardly be recovered from data if only 
time-averaged spectrum were available. In particular, the motion above rnis(a) near a 
rapidly spinning black hole is difficult to distinguish from a descent below rms(O) in the 
non-rotating case. 



3.5 Summary of the new model 



The new ky model is suited for use with the xspec package ( ArnaudI 199d ). Several 



mutations of ky were developed, with the aim to specifically provide different applications, 
and linked with a common ray-tracing subroutine, which therefore does not have to be 
modified when the intrinsic emissivity function in the model is changed. 

When the relativis tic lin e distortions are computed, the new model is more accurate 
than the LAOR model (|l^L991) and faster than KERRSPEC (jMartocchia. Karas Mlttl 



200nh . These are the other two XSPEC models with a sirnilar usage (see also Pariev Sz 



Bromlev 'l998': 'Oierlinski. Maciol ek-Niedzwiecki KbisawalEoOll : iReckwith Donell20oi 



Schnittman &: Bertschineer ■2004lt ) . It is also important to compare the results from fully 



independent relativistic codes since the calculations are sufficiently complex for significant 
differences to arise. Our code is more general than the currently available alternatives in 

XSPEC. 

In our KY model it is possible to choose from various limb darkening/brightening 
laws and thus change the angular distribution of the local emission. The problem of 
directional distribution of the reflected radiation is quite complicated and the matter 
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r\ 1 Vj 1 1 1 It: 






o IJ CC ti dl lillc ii Ulli d Uldt^JY-llUiC U-lol^ 


KYGlline 




non-axisymmetric 




KYFlll 
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non-axisymmetric 


lamp-post fluorescent line 


KYI Irr 




ilUil-clXlby lIllIlcLi IC 


idilip-pUb L v^Uilip Lull 1 clicC LlUll lllUClci, 

physical model of polarization is included 


r\ T n ±1 cl 1 




iioii-axisy mmeLi ic 


v^uiiip Lon 1 eiiecLion wiLii dii incicxeiiL 
powei-iaw loi a uioKeii powei-iawi con- 
tinuum 


KYConv 


convolution 


axisymmetric 


convolution of the relativistic kernel with 


KYClonv 




non-axisymmetric 


intrinsic emissivity across the disc 


KYSpot 


to be used as 


an independent 


time-dependent spectrum of a pattern co- 




code outside xspec 


orbiting with the disc 



' Different model types correspond to the xspec syntax and are defined by the way they act in the 
overall model and form the final spectrum. According to the usual convention in xspec, additive models 
represent individual emission spectral components, which may, for instance, originate in different regions 
of the source. Additive models are simply superposed in the total signal. Multiplicative components (e.g. 
hrefl) multiply the current mode l by an energy-dependent factor. Convolution models modify the model 
in a more non-trivial manner. SeelAinaud (1996) for details. 

Table 3.12: Different versions of the ky model. 



has not been completelv set tled yet (e.g. George fc FabianI 199ll : Zvcki fc Czernvl 19941 : 
Magdziarz &: Zdziarskil 19951 and references cited therein). A specific angular dependence 
is often assumed in models, such as limb darkening /(/Xe) = 1 + 2.06/ic in the LAOR model. 
However, it has been argued that limb brightening may actually occur in the case of 
strong primary irradiation of the disc. This is releva nt for accretion discs near black holes. 
where the effects of emission anisotropy are crucial ( Beckwith fc Don3 20041 : Czernv et al 



20041 ) ■ We find that the spectrum of the inner disc turns out to be very sensitive to the 



adopted angular dependence of the emission, and so the possibility to modify this profile 
and examine the results using ky appears to be rather useful. 

Among its useful features, the ky model allows one to fit various parameters such as 
the black- hole angular momentum (a), the observer inclination angle relative to the disc 
axis (^o)) and the size and shape of the emission area on the disc, which can be non- 
axisymmetric. A straightforward modification of a single subroutine suffices to alter the 
prescription for the disc emissivity, which is specified either by an analytical formula or 
in a tabular form. Our code allows one to change the mesh spacing and resolution for the 
(two-dimensional) polar grid that covers the disc plane, as well as the energy vector (the 
output resolution is eventually determined by the detector in use when the model is folded 
through the instrument response). Hence, there is sufficient control of the (improved) 
accuracy and computational speed. 

Furthermore, KY can be run as a stand-alone program (detached from xspec). In this 
mode there is an o ption for time-variab le sources such as orbiting spots, spiral waves or 
evolving flares (e.g. Czernv et fli.ll20o3 . who applied a similar approach to compute the 
predicted rms variability in a specific flare/spot model). The improved accuracy of the new 
model has been achieved in several ways: (i) photon rays are integrated in Kerr ingoing 
coordinates which follow principal photons, (ii) simultaneous integration of the geodesic 
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deviation equations ensures accurate evaluation of the lensing effect, and (iii) non-uniform 
and rather fine grids have been carefuUy selected. 

Several versions of the routine have been prefabricated for different types of sources 
(Tab. l3?T2|) : (i) an intri nsicallv narrow line produced bv a disc, including the lamp-post 
fluorescent line model (jMartocchia. Karas fc Mattll200oh . (ii) the Compton reflection 
model in two variants - a relativistically blurred Compt on-refiection continuum includin g 
a primary power-law component and a lamp-post model ( Martocchia. Karas Sz Mattlboooh . 
(iii) general convolution models, and (iv) the time-dependent spectrum of an orbiting 
or a free- falling spot. Default parameter values for the line model correspond to those 
in the LAOR model, but numerous options have been added. For example, in the new 
model one is able to set the emission inner radius below the marginally stable orbit, 
Hn < ''^ms(a)- One can also allow a to vary independently, in which case the horizon 
radius, rh(a) = 1 + \/l — a^, has to be, and indeed is, updated at each step of the fit 
procedure. We thus define emission radii in terms of their offset from the horizon. Several 
arguments have been advocated in favour of having rjnT^rms for the disc emission, but 
this possibility has never been tested rigorously against observational data. The set of 
KY-routines introduced above provide the tools to explore black-hole disc models and to 
actually fit for their key parameters, namely, a, Oq, and r[^. 

We have also produced a convolution-type model, kyconv, which can be applied to 
any existing xspec model of intrinsic X-ray emission (naturally, a meaningful combination 
of the models is the responsibility of the user). We remind the reader that kyconv 
is substantially more powerful than the usual convolution models in XSPEC, which are 
defined in terms of a one-dimensional integration over energy bins. Despite the fact that 
KYCONV still uses the standard XSPEC syntax in evaluating the observed spectrum (e.g. 
kyconv(gaussian)), our code performs a more complex operation. It still performs ray- 
tracing across the disc surface so that the intrinsic model contributions are integrated 
from different radii. Thus the kyconv(gaussian) model gives the same results as the 
KYGLINE model if corresponding parameters are set to the same values. The price that 
one has to pay for the enhanced functionality is a higher demand on computational power. 

Other user-defined emissivities can be easily adopted. This can be achieved either by 
using the above mentioned convolution model or by adding a new user-defined model to 
XSPEC. The latter method is more flexible and faster, and hence recommended. In both 
approaches, the ray-tracing routine is linked and used for relativistic blurring. 



Chapter 4 

Applications 



4.1 Seyfert galaxy MCG-6-30-15 

The Seyfert 1 galaxy MCG-6-30-15 is a unique source in which the evidence of a broad and 
skewed Fe K a line has led to a wide acceptance of models with an accreting black hole in 
the n ucleus dXana ka et al}h99^ Iwasawa et al. 1996t iNandra et al. 1997 : Guainazzi et al\ 

m 



19991 : iFabian fc Vaughan.2nn^ Being a nearby AGN (the galaxy redshift is z = 0.0078), 



this source offers an unprecedented opportunity to explore directly the pattern of the 
accretion flow onto the central hole. The Fe K line shape and photon redshifts indicate 
that a large fraction of the emission originates from r < 10 {GM/c^). The mean line 
profile derived from XMM-Newton observations is similar to the one observed previously 
using ASCA. The X-ra y continuum shape in the hard spectral band was well determined 
from BeppoSAX data (jCuainazzi et al.\l99^ . 



To illustrate the new ky model capabilities, we used our code to analyze the mean 
EPIC PN spec trum which we com piled from the long XMM-Newton 2001 campaign (e.g., 
as described in Fabian et a/.ll200^ . The data were cleaned and reduced using standard data 



reduction routines, employing SAS version 5.4.1.^ We summed the EPIC PN data from 
five observations made in the interval 11 July 2001 to 04 August 2001, obtaining a total 
good exposure time of ~ 290 ks (see iFabian alll2002l for details of the observations). 



The energy range was restricted to 3 — 10 keV with 339 energy bins, unless otherwise 
stated, and models were fitted by minimizing the statistic. Statistical errors quoted 
correspond to 90% confidence for one interesting parameter (i.e. Ax^ = 2.706), unless 
otherwise stated. 

No absorber was taken into account; the assumption here is that any curvature in the 
spectrum above 3 keV is not due to absorption, and only due to the Fe K line. We remind 
the reader that the Fe K emission line dominates around the energy 6 — 7 keV, but it has 
been supposed to stretch down to ~ 3 keV or even further. We emphasize that our aim 
here is to test the hypothesis that if all of the curvature is entirely due to the broad Fe K 
feature, is it possible to constrain a of the black hole? Obviously, if the answer to this is 
'no', it will also be negative if some of the spectral curvature between ~ 3 — 5 keV is due 
to processes other than the Fe K line emission. 

We considered Fe Ka and Fe K/3 iron lines (with their rest-frame energy fixed at 
Erest = 6.400 and 7.056 keV, respectively), a narrow-line feature (observed at E'obs ~ 
6.9 keV) modelled as a Gaussian with the central energy, width, and intensity as free 



^See http: //xmin. vilspa. esa. es/sas7| 
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Hn-ni ^b-rh rout-Th Oin "out EW 



1 5.1 ±0.2 - 11.4 ±0.8 - 3.9 ±0.6 258 



2a 
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6 9+°-^ 
o.y_o.4 


q 7+0.3 
^- ' -0.8 
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2b 


0.65 ±0.35 
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4 q+0-4 
^•^-0.3 


941 


2c 


1.23 ±0.06 


4 ± 0.02 




9.2 ±0.2 


3.1 ±0.1 
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Best-fitting values of the important parameters and their statistical errors for models #1-2, described 
in the text. The models include broad Fe Ka and Fe K/3 emission lines, a narrow Gaussian line at 
~ 6.9 keV, and a Compton-reflection continuum from a relativistic disc. These models illustrate 
different assumptions about intrinsic emissivity of the disc (the radial emissivity law does not need to 
be a simple power law, but axial symmetry has been still imposed here). The inclination angle 80 is 
in degrees, relative to the rotation axis; radii are expressed as an offset from the horizon (in GM/(?); 
the equivalent width, EW, is in electron volts. 

1 Two values of the transition radius define the interval (rb-,rb+) where refiection is diminished. 
Table 4.1: Spectral fitting results for MCG-6-30-15 using the ky model. 



parameters, plus a power-law continuum. We used a superposition of two kygline models 
to account for the broad relativistic Ka and K/3 lines, ZGAUSS for the narrow, high-energy 
line (which we find to be centred at 6.86 keV in the source frame, slightly redshifted 
relative to 6.966 keV, the rest-energy of Fe Lya), and kyhIrefl for the Comptonized 
continuum convolved with the relativistic kernel. We assumed that the high-energy narrow 
lin e is non-relativistic. Note that an alternative interpretation of the data (as pointed out 
by Fabian et al. 2002h is that there is not an emission line at 6.9 keV, but He-like 



resonance absorption at ~ 6.7 keV. This would not affect our conclusions. The ratio of 
the intensities of the two relativistic lines was fixed at /k/3/-^Ko = 0.1133 (= 17/150), and 
the iron abunda nce for the Comptq n-refiection continuum was assumed to be three times 
the solar value ( Fabian et al. I2OO2I ). We used the angles ^min = 0° and 0max = 

90° for 

local illumination in KYHlREFL (these values are equivalent to a central illumination of 
an infinite disc). The normalization of the Compton-reflection continuum relative to the 
direct continuum is controlled by the effective reflection 'covering factor', R^, which is the 
ratio of the actual reflection normalization to that expected from the illumination of an 
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Figure 4.1: Confidence contours around the best-fitting parameter values (indicated by a 
cross). Left: the case of a single line-emitting region (model #1 with zero emissivity for 
1^ < '''ms)- Right: the case of a non-monotonic radial emissivity, model i^2h. The joint 
two-parameter contour levels for a versus 6*0 correspond to 68%, 95% and 99% confidence. 



infinite disc. Also rout (outer radius of the disc), Oc (slope of power-law continuum radial 
emissivity in the kyhIrefl model component), as well as Rc, were included among the 
free parameters, but we found them to be only very poorly constrained. Weak constraints 
on Oc and Rc are actually expected in this model for two reasons. Firstly, in these data 
the continuum is indeed rather featureless. Secondly, the model continuum was blurred 
with the relativistic kernel of kyhIrefl, and so the dependence of the final spectrum on 
the exact form of the emissivity distribution over the disc must be quite weak. 

Model #1. Using a model with a plain power-law radial emissivity on the disc, we 
obtained the best-fitting values for the following set of parameters: a of the black hole, 
disc inclination angle 9o, radii of the line-emitting region, < r < rout, the corresponding 
radial emissivity power-law index a of the line emission, as well as the radial extent of 
the continuum-emitting region, its photon index Fc, and the corresponding ac for the 
continuum. Notice that Fc and ac refer to the continuum component before the relativistic 
kernel was applied to deduce the observed spectrum. As a result of the integration across 
the disc, the model weakly constrains these parameters. 

There are two ways to interpret this. A model which is over-parameterized is undesir- 
able from the point of view of deriving unique model parameters from modelling the data. 
However, another interpretation is that a model with a greater number of parameters may 
more faithfully reflect the real physics and it is the actual physical situation which leads 
to degeneracy in the model parameters. The latter implies that some model parameters 
can never be constrained uniquely, regardless of the quality of the data. In practice one 
must apply both interpretations and assess the approach case by case, taking into account 
the quality of the data, and which parameters can be constrained by the data and which 
cannot. If preliminary fitting shows that large changes in a parameter do not affect the fit, 
then that parameter can be fixed at some value obtained by invoking physically reasonable 
arguments pertaining to the situation. 

We performed various fits with the inner edge tied to the marginally stable orbit and 
also fits where rju was allowed to vary independently. Free-fall motion with constant angu- 
lar momentum was assumed below r^s, if the emitting region extended that far. Tab. 14. ll 
gives best-fitting values of the key relativistic line and continuum model parameters for 
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the case in which rjn and a were independent. 

Next, we froze some of the parameters at their best-fit values and examined the 
space by varying the remaining free parameters. That way we constructed joint confidence 
contours in the plane a versus Oq (see left panel of Fig. 14. 1|) . These representative plots 
demonstrate that 6^ appears to be tightly constrained, while a is allowed to vary over a 
large interval around the best fit, extending down to a = 0. 

Model #2. We explored the possibility that the broad-line emission does not conform 
to a unique power-law radial emissivity but that, instead, the line is produced in two 
concentric rings (a 'dual-ring' model). This case can be considered as a toy model for a 
more complex (non-power law) radial dependence of the line emission than the standard 
monotonic decline, which we represent here by allowing for different values of a in the inner 
and outer rings: ain and aout- Effectively, large values of the power-law index represent 
two separate rings. The two regions are matched at the transition radius, r = rb, and 
so this is essentially a broken power law. We explored both cases of continuous and 
discontinuous line emissivity at r = rh- Notice that the double power- law emissivity 
arises naturally in the lamp-post model ( Martocchia. Karas fc Matt] 20001 ) in which the 
disc irradiation and the resulting Fe Ka reflection are substantially anisotropic due to 
fast orbital motion in the inner ring. Although the lamp-post model is very simplified in 
several respects, namely the way in which the primary source is set up on the rotation 
axis, one can expect fairly similar irradiation to arise from more sophisticated schemes of 
coronal flares distributed above the disc plane. Also, in order to provi de a physica l picture 
of the steep emissivity found in XMM-Newton data of MCG-6-30-15, IWilms et al. (2001 i) 
invoked strong magnetic stresses acting in the innermost part of the system, assuming 
that they are able to dissipate a considerable amount of energy in the disc at very small 
radii. Intense self-irradiation of the inner disc may further c ontribute to the effe ct. 

This more complex model is consistent with the findings of lFabian et al. Indeed, 
the fit is improved relative to models with a simple emissivity law because the enormous red 
wing and relatively sharp core of the line are better reproduced thanks to the contribution 
from a highly redshifted inner disc (see Tab. 14. ll model #2a). For the same reason that 
the more complex model reproduces the line core along with the red wing well, the model 
prefers higher values of a and Oo than what we found for case #1. Notice that a — > 1 
implies that all radiation is produced above r^s- Maximum rotation is favoured with both 
a and 9^ appearing to be tightly constrained near their best-fit values. Likewise for the 
continuum radial emissivity indices. There is a certain freedom in the parameter values 
that can be accommodated by this model. By scanning the remaining parameters, we 
checked that the reduced ~ 1 can be achieved also for a going do wn to ^ 0.9 and 



37°. This conclusion is also consistent with the case for large a in lDabrowski et al. 



(19971); however, we actually do not support the claim that the current data require a 
large value of a. As shown below, reasonable assumptions about the intrinsic emissivity 
can fit the data with small a equally well. 

In model #2a, small residuals remain near E ~ 4.8 keV (at about the ~ 1% level), the 
origin of which cannot be easily clarified with the time-averaged data that we employ now. 
The excess is reminiscent of a Doppler horn typical of relativistic line emission from a disc, 
so it may also be due to Fe Ka emission which is locally enhanced on some part of the 
disc. We were able to reproduce the peak by modifying the emissivity at the transition 
radius, where the broken power-law emissivity changes its slope (model #2b). We can 
even allow non-zero emissivity below (the inner ring) with a gap of zero emissivity 
between the outer edge of the inner ring and the inner edge of the outer one. The inner 
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Observed Energy (keV) Observed Energy (keV) 

Figure 4.2: Spectrum and best-fitting model for the XMM-Newton data for MCG-6-30-15 
in which the Fe K Une originates in a dual-ring. The models 2a (left) and 2b (right) are 
shown in comparison. See Tab. l4.ll for the model parameters. The continuum component 
(kyhIrefl) is also plotted. The data points are not unfolded: the spectrum in these 
units was made by multiplying the ratio of measured counts to the counts predicted by 
the best-fitting model and then this ratio was multiplied by the best-fitting model and 
then by E'^. 

ring, rin<r<rb_, contributes to the red tail of the line while the outer ring, rb+<r<rout) 
forms the main body of the broad line. The resulting plot of joint confidence contours of 
9o versus a is shown in Fig. 14.11 fright panel). In order to construct the confidence contours 
we scanned a broad interval of parameters, < a < 1 and < < 45°; here, detail is 
shown only around the minimum region. 

Two examples of the spectral profiles are shown in Fig. 14.21 It can be seen that the 
overall shapes are very similar and the changes concern mainly the red wing of the profile. 
For completeness we also fitted several modifications of the model ^^2 and found that with 
this type of disc emissivity (i.e. one which does not decrease monotonically with radius) 
we can still achieve comparably good fits which have small values of a. The case #2c 
in Tab. 14.11 gives another example. This shows that one cannot draw firm conclusions 
about a (and some of the other model parameters) based on the redshifted part of the line 
using current data. The data used here represent the highest signal-to-noise relativistically 
broadened Fe K line profile yet available for any AGN. 

The differences in between model #1 and models #2a, 2b, and 2c are very large 
(~ 60) for only one additional free parameter, so all variations of model 2 shown are better 
fits than the case #1. We note that formally, model #2a and #2c appear to constrain a 
much more tightly than model #1. The issue here is that the values of a can be completely 
different, with the statistical errors on a not overlapping (compare, for example, model 
7^2a and #2c in Tab. 14. 1|) . This demonstrates that, at least for the parameter a, it is 
not simply only the statistical error which determines whether a small a value is or is not 
allowed by the data. We confirmed this conclusion with more complicated models obtained 
in KY by relaxing the assumption of axial symmetry; additional degrees of freedom do not 
change the previous results. 

One should bear in mind a well-known technical difficulty which is frequently encoun- 
tered while scanning the parameter space of complex models and producing confidence 
contour plots similar to Fig. 14.11 That is, in a rich parameter space the procedure may 
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be caught in a local minimum which produces an acceptable statistical measure of the 
goodness of fit and appears to tightly constrain parameters near the best-fitting values. 
However, manual searching revealed equally acceptable results in rather remote parts of 
the parameter space. Indeed, as the results above show, we were able to find acceptable fits 
with the central black hole rotating either slowly or rapidly ( in terms of the a pa rameter). 
This fact is not in contradiction with previous results (e.g. iFabian et al. 2002^ because 
we assumed a different radial profile of intrinsic emissivity, but it indicates intricacy of 
unambiguous determination of model parameters. We therefore need more observational 
constraints on realistic physical mechanis ms to be able to fit comp licated models to actual 
data with sufficient confide nce |Ballant vne. Ross Sz FabianI 200 ll : iNavakshin & Kazana^ 
20n3:lE6za,nska, et a/.lEooi. 



We have seen that the models described above are able to constrain parameters with 
rather different degrees of uncertainty. It turns out that the more complex type #2 models 
(i.e. those that have radial emissivity profiles which are non-monotonic or even have an 
appreciable contribution fro m r < rms) provide better fits to the data but a physical 
interpretation is not obvious. Ballantvne. Vaughan &: FabianI ( 2003f ) also deduce a dual- 
reflector model from the same data and propose that the outer reflection is due to the disc 
being warped or flared with increasing radius. 
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4.2 Relativistic spectral features from X-ray illuminated 
spots 

Relativistic iron line profiles may provide a powerful tool for measuring the mass of the 
black hole in active galactic nuclei and Galactic black- hole candidates. For this aim, I Stellal 
( 199Qi ) proposed to use temporal changes in the line profile following variations of the il- 
luminating primary source (which at that time was assumed to be l ocated on the disc 
axis for simplicity). Along the same line of thought, iMatt fc Perolal (Il992h proposed to 
employ, instead, variations of the integrated line properties such as equivalent width, cen- 
troid energy and line width. These methods are very similar conceptually to the classical 
reverberation mapping method, widely and successfully applied to optical broad lines in 
AGNs. Sufficiently long monitoring of the continuum and of the line emission is required, 
as well as large enough signal-to-noise ratio. However, the above-mentioned methods have 
not provided many results yet. Even in the best studied case of the Seyfert galaxy MCG-6- 
30-15, the mass estimate is hard to obtain due to the a pparent lack of correlation between 



the line and continuum emission ( Fabian et al. 20021 ). It was also suggested that these 



complica tions are po ssiblv caused by an interplay of complex general relativistic effects 
(Miniutti et al. 20031 ). X-ray spectra from high throughput and high energy resolution 



detectors should resolve the problem of interpretation of observed spectral features. How- 
ever, before such high quality data are available it is desirable to examine existing spectra 
and attempt to constrain physical parameters of the models. 

A simple, direct and potentially robust way to measure the black-hole mass would 
be available if the line emission originates at a given radius and azimuth, as expected 
if the disc illumination is provided by a localized flare just above the disc (possibly due 
to magnetic reconnection), rather than a central illuminator or an extended corona. If 
a resulting 'hot spot' co-rotates with the disc and lives for at least a significant part of 
an orbit, by fitting the light curve and centroid energy of the line flux, the inclination 
angle 9o and the orbit radius could be derived (radius in units of the gravitational radius 
Tg). Further, assuming Keplerian rotation, the orbital period is linked with radius in a 
well-known manner. The equation for the orbital period then contains the black-hole mass 
M, explicitely, and so this parameter can be determined, as discussed later. 

Hot spots in AGN accretion discs were popular for a while, following the finding of 
apparent periodicity in the X-ray emission of the Seyfert 1 galaxy NGC 6814. They, 
however, were largely abandoned when this periodicity was demonstrated to be asso- 



Madeiski et al. 




Iwasawa et all 


1995: 



Periodicities in AGNs were subseq uentlv reported in a few sources 

The fact that they were not 



confirmed in different observations of the same sources is not surprising - quite on the 
contrary, it would be hard to imagine a hot spot surviving for several years. 

Rec ently, the discovery of narrow eni i ssion features in the X-rav spectra of several 



AGNs (jTiiruer et nlhlM] 



Yaqoob et a,l\\20Q^ Turner. Kraemer fc Eeeve'^ 



200l has renewed interest in hot spots. There is a tentative explanation (even if not the 
only one) for these features, typically observed in the 5 — 6 keV energy range, in terms of 
iron emission produced in a small range of radii and distorted by joint action of Doppler 
and gravitational shift of photon energy. Iron lines would be produced by localized flares 
which illuminate the underlying disc surface, producing the line by fluorescence. Indeed, 
the formation of magnetic flares on the disc surface is one of the most promising scenarios 
for the X-ray emission of AGNs. A particularly strong flare, or one with a very large 



54 



CHAPTER 4. APPLICATIONS 



anisotropic emission towards the disc, could give rise to the observed features. Smah 
width of the observed spectral features implies that the emitting region must be small, 
and that it is seen for only a fraction of the entire orbit (either because the flare dies 
out, or because emission goes below detectability, see next section). If the flares co-rotate 
with the disc and if they last for a significant part of the orbit, it may be possible by 
observing their flux and energy variations with phase to determine the orbital parameters, 
and thence M,. 

The basic properties of line emission fro m the innermost regions of an accretion disc 
around a black hole are well-known (see e.g. Fabian et al. 200(1 R.evnolds fc Nowa3 200,4 



for recent reviews). Let us here briefly summarize several formulae most relevant to our 
purposes. 

If r is the orbital radius and a is the dimensionless black-hole angular momentum, the 
orbital period of matter co-rotating along a circular trajectory r = const around the black 
hole is given by (jBardeen. Press fc Teukolskv 1972h 



3 



M, 



To^b = 310 + a) [sec], (4.1) 

as measured by a distant observer. We express lengths in units of the gravitational radius 
rg=GM,/c^=1.48 x lO^^My cm, where My is the mass of the black hole in units of 10^ solar 
ma sses. Angular momentum a (per unit mass) is in geometrized units (0 < a < 1). See 



e.g. lMisner. Thorne fc Wheel er ( 1973.) for useful conversion formulae between geometrized 
and physical units. 

The innermost stable orbit, rms, occurs for an equatorial disc at radius 



1 



r^s = 3 + Z2- [(3-Zi)(3 + Zi + 2Z2)]2, (4.2) 

where Zi = 1 + (1 - a^)^[(l + a)^ + (1 — o)^] and Z2 = (3a^ + Zl)^; r^s spans the range 
of radii from r = 1 (a = 1, i.e. the case of a maximally rotating black hole) to 6 (a = 0, 
a static black hole). Rotation of a black hole is believed to be limited by a n equilibrium 
value a=0.998 because of the capture of photons from the disc (jThornelll974h . This would 
imply rms=1.23. Different specific models of accretion can result in somewhat different 
limiting values of a and the corresponding r^^. Notice that in the static case, the radial 
dependence Torb(?^)|a=o is identical to that in purely Newtonian gravity. 

In order to compute a synthetic profile of an observed spectral line one has to link 
the points of emission in the disc with corresponding pixels in the detector plane at spa- 
tial infinity. This can be achieved by solving the ray-tracing problem in curved space- 
time of the black hole. Appropriate methods were discussed by several authors; see 
lR,evn olds & Nowak (20q3) for a recent review and for further references. This way one 
finds the redshift factor, which determines the energy shift of photons, the lensing effect 
(i.e. the change of solid angle due to strong gravity), and the effect of aberration (which 
influences the emission direction of photons from the disc; this must be taken into account 
if the intrinsic emissivity is non-isotropic). We con sider these effects in our computa tions, 
assuming a rotating (Kerr) black- hole space-time ( Misner. Thorne &: Wheeler 19731 ). We 



also consider the time of arrival of photons originating at different regions of the disc plane. 
Variable travel time results in mutual time delay between different photons, which can be 
ignored when analyzing time-averaged data but it may be important for time-resolved 
data. 
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Figure 4.3: Line flux and centroid energy as functions of the orbital phase of a spot, for 
three values of angular momentum (a = 1, 0.9, and 0) and three inclination angles (^o = 
30°, 60°, 85°). The centre of the spot is located at radial distance r, which corresponds to 
the last stable orbit rms(a) for that angular momentum plus a small displacement given by 
the spot radius, dr. The intrinsic energy of the line emission is assumed to be at 6.4 keV 
(indicated by a dotted line). Prograde rotation is assumed. Time is expressed in orbital 
periods. 
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Figure 4.4: The same as the previous figure, but with a = and r = 10 (top), 20 (middle), 
and 50 (bottom). 
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Figure 4.5: Line profiles integrated over twelve consecutive temporal intervals of equal 
duration. Each interval covers 1/12 of the orbital period at corresponding radius. As 
explained in the text, top-left frame of each panel corresponds to the spot being observed 
at the moment of passing through lower conjunction. Energy is on abscissa (in keV). 
Observed photon flux is on ordinate (arbitrary units, scaled to the maximum flux which 
is reached during the complete revolution of the spot). Notice the occurrences of narrow 
and prominent peaks which appear for relatively brief fraction of the total period. 
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Figure 4.6: The same as in previous figure, but with r = 10 (top), 20 (middle), and 50 
(bottom). The black hole was assumed non-rotating, a = in this figure. 
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Assuming purely azimuthal Keplerian motion of a spot, one obtains for its orbital 
velocity (with respect to a locally non-rotating observer at corresponding radius r): 

= (4.3) 



A (r3/2 + a) 

In order to derive time and frequency as measured by a distant observer, one needs to 
take into account the Lorentz factor associated with this orbital motion, 

(r3/2 + a) 



rV4 v^r3/2 _ 3^1/2 _^ 2a Vr^ + a?r + 2a^ 

The corresponding angular velocity of orbital motion is 17 = (r3/2 + a)~^, which also 
determines the orbital period in eq. H4.1|) . The redshift factor g and the emission angle •& 
(with respect to the normal direction to the disc) are then given by 

= arccos (4.5) 



" e-r-3/2^' 

where i3 = 1 + ar^3/2^ — y — 3r~^ + 2ar^3/2. ^ ^ ^j,g constants of motion connected 
with the photon ray in an axially symmetric and stationary space-time. 

For practical purposes formula 1)4. 1|) with a = is also accurate enough in the case of a 
spinning black hole, provided that r is not very small. For instance, even for r = 6 (the last 
stable orbit in Schwarzschild metric), Toxhii^ms) calculated for a static and for a maximally 
rotating (a = 1) black hole differ by about 6.8%. The relative difference decreases, roughly 
linearly, down to 1.1% at r = 20. This implies that eq. (|4.1() can be used in most cases to 
estimate the black-hole mass even if the angular momentum is not known (deviations are 
relevant only for r < 6, when the radius itself can be used to constrain the allowed range 
of a). 

Various pseudo-Newtonian formulae have been devised for accreting black holes to 
model their observatio nal properties, which are c onnected with the orbital motion of sur- 



roundin g matter (e.g. Abramowicz et a,l. 199fil : Artemova. Biornsson fc Novikov 199fil : 
ISemerak fc Karaill99g ). Although this approach is often used and found to be practical, 



we do not employ it here because error estimates are not possible within the pseudo- 
Newtonian scheme. 

Due to Doppler and gravitational energy shift the line shape changes along the orbit. 
Centroid energy is redshifted with respect to the rest energy of the line emission for most 
of the orbit. Furthermore, light aberration and bending cause the flux to be strongly 
phase-dependent. These effects are shown in Figs. I4.3H4.4I In these plots, the arrival time 
of photons is defined in orbital periods, i.e. scaled with To^Y,{r;a). The orbital phase of 
the spot is of course linked with the azimuthal angle in the disc, but the relation is made 
complex by time delays which cannot be neglected, given the large velocities of the orbiting 
matter and frame-dragging effects near the black hole. Here, zero time corresponds to the 
moment when the centre of the spot was at the nearest point on its orbit with respect to the 
observer (a lower conjunction). The plots in Fig. 14. HI refer to the case of a spot circulating 
at the innermost stable orbit rjns(a) for a = 0, 0.9 and 1. The effect of black-hole rotation 
becomes prominent for almost extreme values of a; one can check, for example, that the 
difference between cases a = and a = 0.5 is very small. 

Worth remarking is a large difference in the orbital phase of maximum emission be- 
tween the extreme case, a — > 1, in contrast to the non-rotating case, a ^ 0. The reason 
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a=0, 1=1^^ a=0.9, i^i^^ a=1 , i=i^ 




Figure 4.7: Time-averaged synthetic spectra in terms of photon flux (in arbitrary units) 
versus energy (in keV). These profiles represent the mean, background-subtracted spectra 
of the Fe Ka iron-line originating from spots at different radii. Top panels correspond 
to r = Tins and a = (left), a = 0.9 (middle), and a = 1 (right). In bottom panels we 
fix a = and choose r = 10, 20, and 50, respectively (other values of a give very similar 
profiles). Three consecutively increasing values of observer inclination Oq are shown, as 
indicated inside the frames. 



is that for large a the time delay and the effect of frame-dragging on photons emitted 
behind the black hole are very substantial. It is also interesting to note that, for very 
high inclination angles, most of the flu x comes from the f a r side of the disc, due to ver\ 



strong light bending, as pointed out bv lMatt et al. I (|l992h : iMatt. Perola fc St!^ 
and examined further by many authors who performed detailed ray-tracing, necessary to 
determine the expected variatio ns of the l ine flux and shape. A relatively simple fltting 
formula has also been derived (jKaraslll996h and can be useful for practical computations. 

Three more orbits (centred at r = 10, 20 and 50) are shown for a = (Fig. 14.411 . As 
said above, at these radii differences between spinning and static black holes are small. 
Indeed, it can be verified that the dependence on a is only marginal if r > 20, and so it 
can be largely neglected for present-day measurements. 

In Figs. I4.5H4.6I we show the actual form of the line profiles for the same sets of 
parameters as those explored in Figs. I4.3H4.4I The entire revolution was split into twelve 
different phase intervals. The intrinsic flux / is assumed to decrease exponentially with 
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the distance dr from the centre of the spot (i.e. log/ cx — [«;dr/r]^, where r is the location 
of the spot centre and k ~ 10 is a constant). The illumination is supposed to cease 
at dr = 0.2r, which also defines the illuminated area in the disc. Let us remark that 
we concentrate on a spectral line which is intrinsically narrow and unresolved in the rest 
frame of the emitting medium. Such a line can be produced by a spot which originates due 
to s harply localized illumination by flares, as proposed and discussed by various authors 
(e.g. Haardt. Maraschi fc GhiselliriilllQQ^ Poutanen FabianI 19991 : Merloni &: FabianI 
Very recently. ICzerny et al\ (|2004l ) have examined the induced rms variability in 



the flare/spot model with relativistic effects. In many cases, and especially for small radii 
and intermediate to large inclination angles, the line emission comes from a relatively 
minor fraction of the orbit. This implies in practice that for observations with moderate 
signal-to-noise ratios, only a narrow blue horn can be visible, and only for a small part of 
the orbit. These large and rapid changes of the line shape get averaged when integrating 
over the entire orbit, and so an important piece of information is missing in the mean 
spectra. The line profiles integrated over the whole revolution are shown in Fig. 14.71 
Effectively, the mean profile of a spot is identical to the profile of an annulus whose radius 
is equal to the distance of the spot centre and the width is equal to the spot size. 

We discussed the possibility that the narrow features in the 5 — 6 keV range, recently 
discovered in a few AGNs and usually interpreted as redshifted iron lines, could be due to 
illumination by localized orbiting spots just above the accretion disc. If this is indeed the 
case, these features may provide a powerful and direct way to measure the black-hole mass 
in active galactic nuclei. To achieve this aim, it is necessary to follow the line emission 
along the orbit. The orbital radius (in units of rg) and the disc inclination can be inferred 
from the variations of the line flux and centroid energy. Furthermore, M, can be estimated 
by comparing the measured orbital period with the value expected for the derived radius. 
We must point out, however, that present-day X-ray instruments do not have enough 
collecting area to perform this task accurately. This capability should be achieved by the 
planned high-performance X-ray missions such as Constellation- X and Xeus. 



62 



CHAPTER 4. APPLICATIONS 



4.3 Polarization signatures of strong gravity in AGN accre- 
tion discs 



Accretion discs in central regions of active galactic nuclei are subject to strong external 
illumination originating from some kind of corona and giving rise to specific spectral fea- 
tures in the X-ray band. In particular, the K-shell lines of iron are found to be prominent 
around 6 — 7 keV. It has been shown that the shape of the intrinsic spectra must be further 
modified by the strong gravitational field of the central mass, and so X-ray spectroscopy 
could allow us to explore the innermost r egions o f accretion flo ws near supermassive black 
holes (for recent reviews see lFabian et al.. 2000: Revnolds &: Nowak.2003. ). Similar mech- 
anisms operate also in some Galactic black-hole candidates. 

A rather surprising result from recent XMM-New ton observations is t hat relativistic 
iron lines are not as common as previously believed, see lBianchi et al. I (l20n4h and references 
therein and also lYaaoob et al, (,2003. ) ■ This does not necessarily mean that the iron line 
is not produced in the innermost regions of accretion discs. The situation is likely to 
be more complex than in simple, steady scenarios, and indeed some evidence for line 
emission arising from orbiting spots is present in the time-resolved spectra of a few AGNs 
(jPovciak et a/.. .2004ai ) . Even when clearly observed, relativistic lines behave differently 
than expected. The best example is the puzzling lack o f correlation between line and 
continuum emission in MCG~6-30-15 dFabian et fli.ll2002h . unexpected because the very 
broad line profile clearly indicates that the line originates in the innermost r e gions of the 
accretion disc, hence very close to the illuminating source. iMiniutti etld\ (|2003h have 
proposed a solution to this problem in terms of an illuminating source moving along the 
black-hole rotation axis or very close to it. 

In this section we show that polarimetric studies could provide additional information 
about accretion discs in a strong gravity regime, which may be essential to discrimi- 
nate between different possible geometries of the source. The idea of using polarimetry 
to gain additional information about accre ting compact objects is not a new one. In 
this context it was proposed bv lReed ( 19751 ) that polarized X-rays are of high relevance. 
Pozdnvakov. Sobol' fc SunvaevI (|l97iy studied spectral profiles of iron X-ray lines that re- 
sult from multiple Compton scattering. Later on, various influences affecting polarization 
(due to magnetic fie lds, absorp tion as well as strong gravity) were examined for black- 
hole accretion discs ( Aeoll 199ffl . Temporal variations of p olarization were also di scussed. 



in^articular the cas^TbitTng spots near a black hole JCouuors. Pira. X. Sta.kd98r^: 



Bao Wiita & H adravJ l99fi). With the promise of new polarimetric detectors ([Costa et al 
200lh . quantitative examination of specific models becomes timely. 



Since the reflecting medium has a disc-like geometry, a substantial amount of linear 
polarization is expected in the resulting spectrum because of Compton scattering. Po- 
larization properties of the disc emission are modified by the photon propagation in a 
gravitational field, providing additional information on its structure. Here we calculate 
the observed polarization of the reflected ra diation assuming t he lamp -post model for 
the s tationary power-law illuminating source ( Martocchia Sz Matti.l99& : .Petrucci Sz Henri 
1223) • We assume a rotating (Kerr) black hole as the only source of the gravitational 
field, having a common symmetry axis with an accretion disc. The disc is also assumed 
to be stationary and we restrict ourselves to the time-averaged analysis. In other words, 
we examine processes that vary at a much slower pace than the light-crossing time at 
the corresponding radius. Intrinsic polarization of the emerging light can be computed 
locally, assuming a plane-parallel scattering layer which is illuminated by light radiated 



4.3. Polarization signatures of strong gravity in AGN accretion discs 



63 



=30°, rin=6, rout=400 



=60°, r|n=6, rout='*00 





Height 




2 




6 

15 




100 






'/'ll 








4 6 8 10 12 14 
Energy [keV] 

9,=30°, r;„=6, r.„,=400 



Height refiected+direct 

2 

6 

15 

100 





Height 
2 




6 




15 




100 




^ _ 










4 6 


8 10 12 14 
Energy [l<eV] 


^° 


-Rn° r -R r -400 




rellected Height 

2 




6 




15 




100 










4 6 


8 10 12 14 
Energy [l<eV] 




=60°, r|n=6, rout=400 


Height 

2 


refiected+direct 


6 




15 




100 - 





60=80°, r|n=6, rout=400 

^ ' ' Height 

2 

6 

15 

100 

6 8 10 12 14 
Energy [keV] 

8o=80", ri„=6, r„„,=400 

refiected Height 

2 

6 

15 
100 



6 8 10 12 
Energy [keV] 

e„=80°, rin=6, r„,,i=400 



Height refiected+direct 

2 



15 
100 



4 6 8 10 12 14 
Energy [keV] 



6 8 10 12 14 
Energy [keV] 



8 10 12 14 

Energy [keV] 



Figure 4.8: Energy dependence of polarization angle (top panels) and polarization degree 
(middle panels) due to reflected radiation for different observer's inclination angles (^o = 
30°, 60° and 80°) and for different heights of the primary source {h = 2, 6, 15 and 100). 
Polarization degree for reflected plus direct radiation is also plotted (bottom panels). The 
emission comes from a disc within = 6 and rout = 400. Isotropic primary radiation 
with photon index F = 2 and angular momentum of the central black hole a = 0.9987 
were assumed. 




from the primar y source. This problem w as studied ex tensively in various approximations 
(e.g. ICha ndras ekhail 19601: Sunva cv & Tr 
computations (jMatt. Perola fc Pirdlll99 



Here we employ the Monte Carlo 
and thus we find the intrinsic emis- 
sivity of an illuminated disc. Then we integrate contributions to the total signal across the 
disc emitting region using a general relativistic ray-tracing technique described in previ- 
ous chapters and we compute the polarization angle and degree as measured by a distant 
observer (see Section [2.31 and equations therein). We show the polarization properties of 
scattered light as a function of model parameters, namely, the height z = h of the primary 
source on the symmetry axis, the dimensionless angular momentum a of the black hole, 
and the viewing angle 9o of the observer. 

In the first set of figures (Figs. I4.8H4.9() we show the energy dependence of polarization 
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Figure 4.9: Same as in the previous figure but for disc starting at Hn = 1.20 . 



angle and degree due to reflected and reflected plus direct radiation for different inclination 
angles and different heights of the primary source. One can see that the polarization of 
reflected radiation can be as high as thirty percent or even more for small inclinations 
and small heights. Polarization of the reflected radiation does not depend on energy very 
much except for the region close to the iron edge at approximately 7.2 keV, where it either 
decreases for small inclinations or increases for large ones. 

In order to compute observable characteristics one has to combine the primary power- 
law continuum with the reflected component. The polarization degree of the resulting 
signal depends on the mutual proportion of the two components and also on the energy 
range of an observation. The overall degree of polarization increases with energy (see 
bottom panels in Figs. l4.8H4.9j) due to the fact that the intensity of radiation from the 
primary source decreases exponentially, the intensity of the reflected radiation increases 
with energy (in the energy range 3—15 keV) and the polarization of the reflected light alone 
is more or less constant. In our computations we assumed that the irradiating source emits 
isotropically and its light is affected only by gravitational redshift and lensing, according to 
the source location at z = /i on axis. This results in a dilution of primary light by factor 
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20 40 60 80 100 20 40 60 80 100 

Height Height 

Figure 4.10: Polarization degree and angle due to reflected radiation integrated over the 
whole surface of the disc and propagated to the point of observation. Dependence on 
height h is plotted. Left panel: = 6; right panel: rm = 1.20. In both the panels the 
energy range was assumed 9 — 12 keV, the photon index of incident radiation F = 2, the 
angular momentum a = 0.9987. 




Figure 4.11: Polarization degree and angle as functions of fio (cosine of observer inclination, 
/io = corresponds to the edge-on view of the disc). The same model is shown as in the 
previous figure. 
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Figure 4.12: Net polarization degree of the total (primary plus reflected) signal as a 
function of h. Left panel: Tin = 6; right panel: r-m = 1-20. The curves are parametrized 
by the corresponding energy range. 




Figure 4.13: Net polarization degree of the total (primary plus reflected) signal as a 
function of ^o- The same model is shown as in the previous figure. 
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^g\{h,9o)l]x{h,9o), where (7h = — 2h/ {a? + /i^) is the redshift of primary photons 
reaching directly the observer, is the corresponding lensing factor. Here, the redshift is 
the dominant relativistic term, while lensing of primary photons is a few percent at most 
and it can be safely ignored. Anisotropy of primary radiation may further attenuate or 
amplify the polarization degree of the final signal, while the polarization angle is rather 
independent of this influence as long as the primary light is itself unpolarized. 

The polarization of scattered light is also shown in Fig. I4.1U1 where we plot the polar- 
ization degree and the change of the polarization angle as functions of h. Notice that in 
the Newtonian case only polarization angles of 0° or 90° would be expected for reasons 
of symmetry. The change in angle is due to gravitation for which we assumed a rapidly 
rotating black hole. The two panels in the figure correspond to different locations of the 
inner disc edge: = 6 and nn = 1.20, respectively. The curves are strongly sensitive to 
r\n and /i, while the dependence on rout is weak for a large disc (here rout = 400). Sensitiv- 
ity to Tin is particularly appealing if one remembers the practical difficulties in estimating 
Hn by fitting spectra. The effect is clearly visible even for h ^ 20. Graphs corresponding 
to rin = 6 and a = 0.9987, resemble, in essence quite closely, the non-rotating case (o = 0) 
because dragging effects are most prominent near the horizon. 

Fig. 14.111 shows the polarization degree and angle as functions of the observer's incli- 
nation. Again, by comparing the two cases of different one can clearly recognize that 
the polarization is sensitive to details of the flow near the inner disc boundary. 

The dependence of the polarization degree of overall radiation (primary plus reflected) 
on the height of the primary source and the observer inclination in different energy ranges 
is shown in Figs. I4.12H4.151 

In this section we examined the polarimetric properties of X-ray illuminated accretion 
discs in the lamp-post model. From the figures shown it is clear that observed values of 
polarization angle and degree are rather sensitive to the model parameters. The approach 
adopted provides additional information with respect to traditional spectroscopy and so 
it has great potential for discriminating between different models. It offers an improved 
way of measuring rotation of the black hole because the radiation properties of the inner 
disc region most likely reflect the value of the black-hole angular momentum. 

While our calculations have been performed assuming a stationary situation, in reality 
it is likely that the height of th e illuminating s ource changes with time, and indeed such 
variations have been invoked bv lMiniutti et al. to explain the primary and reflected 

variability patterns of MCG-6-30-15. A complete time-resolved analysis (including all 
consequences of the light travel time in curved space-time) is beyond the scope of this 
section and we defer it to future work, assuming that the primary source varies on a 
time-scale longer than light-crossing time in the system. This is also a well-substantiated 
assumption from a practical point of view, since feasible techniques will anyway require 
sufficient integration time (i.e. order of several ksec). Once full temporal resolution is 
possible, the analysis described above can be readily extended. Here, it suffices to note 
that a variation of h implies a variation of the observed polarization angle of the reflected 
radiation. As it is hard to imagine a physical and/or geometrical effect giving rise to the 
same effect, time variability of the polarization angle can be considered (independently of 
the details) a very strong signature of strong-f ield general relati vitv effects at work. 

New generation photoelectric polarimeters ( Costa et a/.ll200lh in the focal plane of large 
area optics (such as those foreseen for Xeus) can probe the polarization degree of the order 
of one percent in bright AGNs, making polarimetry, along with timing and spectroscopy, 
a tool for exploring the properties of the accretion flows in the vicinity of black holes. 



Summary and future prospects 



In this thesis we have developed a computational tool for modelling spectral features of 

X-ray sources in a strong gravitational field. 

The following six transfer functions have been computed for light rays emitted from 
the equatorial plane of the Kerr black hole and received by an observer at infinity: the 
gravitational and Doppler shift, lensing, emission angle, relative time delay, change of the 
polarization angle and azimuthal emission angle. The values of these functions for different 
angular momenta of the black hole and inclination angles of the observer in Kerr ingoing 
coordinates have been stored in the form of tables in a FITS file. In computations, where 
parameters of the motion of the matter were needed, we have assumed a co-rotating Kep- 
lerian disc above the marginally stable orbit and freely falling matter below it. Graphical 
representations of the results have been made in the form of an atlas of contour figures of 
these functions. 

For the modelling of the spectra of an accretion disc, several general relativistic codes 
have been developed. Sonic of them can be used inside a more general spectral fitting 
package XSPEC for routine fitting of the data measured by X-ray satellites. These are the 
models for the relativistic emission lines (kygline, kygIline and kyfIll), Compton 
reflection (kylIcr and kyhIrefl) and for the general use as a relativistic extension of 
existing models (the convolution models kyconv and KYClONv). Other components of 
the tool have extended features for stTidying the non-stationary emission from accretion 
discs (kygspot) and polarimetry (kylIcr). 

The newly developed codes have been employed in several applications. Firstly, we 
used our new KY models inside XSPEC and fitted the data from XMM-Newton observations 
of Seyfert galaxy MCG-6-30-15. Then we simulated emission from an X-ray illuminated 
spot orbiting near a black hole as an application on non-stationary emission from the 
disc. And finally, the polarization from an accretion disc illuminated by a primary source 
located above the black hole was computed. 

There are several possibilities of how one can proceed further. We would like to: 

• enhance the code so that it would be possible to investigate the self-irradiation of 
the disc, 

• add more models to the ky package, first of all a model for the thermal emission, 

• expand the code so that non-Keplerian discs can be studied as well, 

• perform more polarization studies with the non-stationary lamp-post model, 

• perform timing analysis of the emission from an orbiting spot in different energy 
bands, 

• fit more data of X-ray sources acquired by X-ray satellites. 

Further it would be possible to investigate applications of our code to quasi-periodic os- 
cillations. Although our code is a two-dimensional one it might be interesting to compute 
how oscillations of the disc in the equatorial plane can affect the observed spectra at 
infinity. 



Appendix A 

Kerr space-time 



In this Appendix as well as everywhere else in this thesis we use units where GM, = c = 1 
(M, is the mass of the central black hole) and we assume that the angular momentum a 
of the black hole is positive. 



A.l Kerr metric 

The Kerr metric in the Boyer-Lindquist coordinates {t,r,0,(p) is 
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where = + cos^O, A = — 2r + and A = {r"^ + a?Y ~ sm^O. 

Let's define special Kerr ingoing coordinates (t, u, fi, (p) by the following tetrad vectors 
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The advantage of these coordinates, where u = and jl = cos 6, is that the Kerr 
metric is not singular on the horizon of the black hole and spatial infinity (r — > oo) is 
brought to a finite value (u ^ 0). Another advantage for numerical computations is that 
we get rid of the cosine function. 

The relationship between the Boycr-Lindquist coordinate ip and the Kerr ingoing co- 
ordinate ipK, for geodesies coming from infinity to the equatorial plane, can be expressed 
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in the following way: 
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with r± = li: \/l — being inner (— ) and outer (+) horizon of the black hole. 

Another useful way to express the transformation between the Boyer-Lindquist and 
special Kerr ingoing coordinates is by matrices of transformation 
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The Kerr metric in special Kerr ingoing coordinates is 



\ 



5^ = p^/r^ 
,2 -2 An 



-(l-f) 


1 







1 








a(l-A2) 















a{l-/i2) 





^(i-A^) 


■u^/i^, A = 


A/r2 = 


1 - 2n + 


v? and ^ = 



(A.7) 



,2 ~2\2 



a2n2A(l -/i2). 

A. 2 Light rays in Kerr space-time 

The four-momentum = of photons t ra velling in Kerr space-time in the Boyer 



Lindquist coordinates is (see e.g. ICarteiiri968l or Misner. Thorne &: Wheeler 19731 ) 
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A. 3. Keplerian disc in Kerr space-time 
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where I = a(l— /Xq)^/^ = asin^o and q"^ = P'^+fi^{a'^ — a^) are Carter's constants of motion 
with a and P being impact parameters measured perpendicular and parallel, respectively, 
to the spin axis of the black hole projected onto the observer's sky. Here we define a to 
be positive when a photon travels in the direction of the four-vector ^ at infinity, and 

P to be positive if it travels in the direction of — ^ at infinity. The parameter (and 
/io = COS ^o) defines a point at infinity through which the light ray passes (we consider only 
light rays coming to or from the observer at infinity). Furthermore, we have denoted the 
sign of the radial component of the momentum by Rsgn and the sign of the ^-component 
of the momentum by Gggn- We have chosen an affine parameter A' along light geodesies 
in such a way that the conserved energy is normalized to —pt = 1. 

The four-momentum transformed into special Kerr ingoing coordinates is 
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l + 2u-{P + q^)u^ 



u{2 — aul) + UsgnVU 



(A.12) 

(A.13) 
(A.14) 
(A.15) 



with U = l + {a'^-P-q'^)u'^ + 2[{a-lf + q'^]u^-a\'^u^ and M = q'^ + {a^-P-q^)p^- 
The sign of the u- and /i-components of the four- momentum is denoted by C/sgn and Msgn, 
respectively. Note that in these coordinates the four-momentum is not singular on the 
horizon as opposed to the expressions in the Boyer-Lindquist coordinates. 



A. 3 Keplerian disc in Kerr space-time 

Matter m oves a long free stable circular orbits with a rotational velocity (see e.g. Novikov 
k ThorneEszi) 

dt r-*/^ -I- a 

in a Keplerian disc. The disc resides in the equatorial plane of the black hole. Individual 
components of the four-velocity of the Keplerian disc in the Boyer-Lindquist coordinates 
are 

= + , (A.17) 

^\/ — 3r + 2ay/r 

V = 0, (A.18) 
f/^ = 0, (A.19) 

U'f' = ^ (A.20) 

y/ j.(j.2 _ _|_ 2ay/r) 

There is no free stable circular orbit below the marginally stable orbit defined by 



3 + Za - V(3-Zi)(3 + Zi + 2Z2)] (A.21) 
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with = 1 + (1 - a2)V3[(i + a)i/3 + (i _ a)i/3] and Z2 = ^/3a?l^. We suppose that 
below this orbit the matter is in a free fall down to the horizon. Thus the matter conserves 
its specific energy —Ut and its specific angular momentum 

-Ut{r<r^) ^ ^ -Utir^s) = ^^"^ + ^ (a.22) 

-2 I -2 2a,/nr 



When we also consider the normalization condition for the four- velocity, U^Ufj, = —1, 
then we get the following expressions for its contravariant components 

U\r<r^,) = -L{[r(r2 + a2) + 2a2]S^, -2aL„,J, (A.24) 
1 



t/"(r<r^s) = ^[r(r2 + a^) + 2a^]El, - ^aE^,L^, - (r - 2)LL - rA , (A.25) 

C/^(r<r^s) = 0, (A.26) 

i7^(r<r^s) = ^[2a^ms + (r-2)L„J. (A.27) 

In our calculations we use the following local orthonormal tetrad connected with the 
matter in the disc 

e(i)^ = f/^, (A.28) 

e(,)^ ^ ^==^=== [ [/'^(C/t, 0, C/^) + (0, 1, 0, 0) ] , (A.29) 
y A(i + U Ur) 

6(5)^ = (0, 0, r, 0), (A.30) 



{-V^, 0, 0, [/*) . (A.31) 



1 + U^Ur 

A. 4 Initial conditions of integration 

In this part of the Appendix we will sum up the initial conditions for integration of the 
equation of the geodesic 

^ = %+^^-rP''P^, fi = i,u,fi,^ (A.32) 
and equation of the geodesic deviation 
d^Y;^ dK^ 

^ + 2r^^K^ + P^,^^pV>^^ = 0, ^, = t,u,pL,^, j = l, 2 (A.33) 

which we solve in order to compute transfer functions over the accretion disc (see Chapter^ 
for details). Here = is a four- momentum of light, A' is an affine parameter for which 
the conserved energy —pt = 1, Y^^ are two vectors characterizing the distance between 
nearby geodesies and F^^ are Christoffel symbols for the Kerr metric. We solve these 
equations in special Kerr ingoing coordinates t, u, fi, ip and we start integrating from the 
following initial point: 

ti = , u\ = 10~^^ , /ii = /io = cos 6*0 , = , (A. 34) 



A. 5. Christoffel symbols 
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where Oq is the inclination angle of the observer at infinity. The initial values of the four- 
momentum of light were defined by this point and eqs. (|A.12p - (|A.15|) with Uggn = +1, 
Msgn = for /3 = and Afsgn = sign(/3) for /? / 0. 
Four-vectors Y"/* and Y2 have the initial values 



1 ' 



and their derivatives have the initial values {dY^/dX' = dY^ / dx'^) 



J ' ^ i 

dy* Hy* 

d# = d# = 

HV" dy" 

dA' " dA' " 

dY.^. „ / dY^, „ 

dy^f dy,f 



(A.35) 



(A.36) 



The initial vectors Yy^ and l^'j' were chosen in such a way that in the initial point, where 
we suppose the metric to be a flat-space Minkowski one, they are perpendicular to each 
other and to the four-momentum of light, are space-like and have unit length. We kept 
only the largest terms in u\ in eqs. ()A.35|1 and ()A.36|) . 



A. 5 Christoffel symbols 

Here we show the Christoffel symbols and their first derivatives for the Kerr metric in 
special Kerr ingoing coordinates. They are needed for the numerical integration of the 
equation of the geodesic and geodesic deviation. 
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au{l - ji^)[a^u^{l - - a^i?t?) - ~p'^\/~p^ 
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r^- . =0 

UU,fi " 
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r^^,^ = 2uAA[^2 - a2u2(i - a2)]/[(i - A') V] 
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Description of FITS files 



B.l Transfer functions in KBHt able sNN. fits 

The transfer functions are stored in the file KBHtablesNN . f its as binary extensions and 
parametrized by the value of the observer inclination angle Oq and the horizon of the 
black hole r^- We found parametrization by rh more convenient than using the rotational 
parameter a, although the latter choice may be more common. Each extension provides 
values of a particular transfer function for different radii, which are given in terms of r — rh, 
and for the Kerr ingoing axial coordinates </?k- Values of the horizon rh, inclination 9o, 
radius r — rh and angle ipK, at which the functions are evaluated, are defined as vectors 
at the beginning of the FITS file. 

The definition of the file KBHtablesNN. fits: 

0. All of the extensions defined below are binary. 

1. The first extension contains six integers defining which of the functions is present 
in the tables. The integers correspond to the delay, ^-factor, cosine of the local 
emission angle, Icnsing, change of the polarization angle and azimuthal emission 
angle, respectively. Value means that the function is not present in the tables, 
value 1 means it is. 

2. The second extension contains a vector of the horizon values in GM/c^ (1.00 < rh < 
2.00). 

3. The third extension contains a vector of the values of the observer's inclination angle 
Oo in degrees (0° < 6*0 < 90°, 0° - axis, 90° - equatorial plane). 

4. The fourth extension contains a vector of the values of the radius relative to the 
horizon r — rh in GM/ . 

5. The fifth extension contains a vector of the values of the azimuthal angle (pK in 
radians (0 < ipK < Svr). Note that ipK is a Kerr ingoing axial coordinate, not the 
Boyer-Lindquist one! 

6. All the previous vectors have to have values sorted in an increasing order. 

7. In the following extensions the transfer functions are defined, each extension is for 
a particular value of rh and 9o- The values of rh and 60 are changing with each 
extension in the following order: 
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rh[l] X 0o[2], 
rh[l] X 0,[3], 



rh[2] X 0o[l], 
rh[2] X eo[2], 
rh[2] X 0o[3], 

8. Each of these extensions has the same number of columns (up to six). In each 
column, a particular transfer function is stored - the delay, (7-factor, cosine of the 
local emission angle, lensing, change of the polarization angle and azimuthal emission 
angle, respectively. The order of the functions is important but some of the functions 
may be missing as defined in the first extension (see 1. above). The functions are: 

delay - the Boyer-Lindquist time in GM/c" that elapses between the emission of 
a photon from the disc and absorption of the photon by the observer's eye at 
infinity plus a constant, 

^-factor - the ratio of the energy of a photon received by the observer at infinity to 
the local energy of the same photon when emitted from an accretion disc, 

cosine of the emission angle - the cosine of the local emission angle between the 
emitted light ray and local disc normal, 

lensing - the ratio of the area at infinity perpendicular to the light rays through 
which photons come to the proper area at the disc perpendicular to the light 
rays and corresponding to the same flux tube, 

change of the polarization angle in radians - if the light emitted from the disc is 
linearly polarized then the direction of polarization will be changed by this 
angle at infinity - counter-clockwise if positive, clockwise if negative (we are 
looking towards the coming emitted beam); on the disc we measure the angle 
of polarization with respect to the "up" direction perpendicular to the disc 
with respect to the local rest frame; at infinity we also measure the angle of 
polarization with respect to the "up" direction perpendicular to the disc - the 
change of polarization angle is the difference between these two angles, 

azimuthal emission angle in radians - the angle between the projection of the three- 
momentum of an emitted photon into the disc (in the local rest frame co- moving 
with the disc) and the radial tetrad vector. 

For mathematical formulae defining the functions see eqs. (|1.5j) -- (|1.7|) . (|1.1())1 - (|1.12|) 
and (fTTH)) in Chapter U] 

9. Each row corresponds to a particular value of r — rh (see 4. above). 

10. Each element corresponding to a particular column and row is a vector. Each element 
of this vector corresponds to a particular value of (/?k (see 5. above). 

We have pre-calculated three sets of tables - KBHtablesOO . f its, KBHtables50 . f its 
and KBHtables99 . f its. All of these tables were computed for an accretion disc near a 
Kerr black hole with no disc corona present. Therefore, ray-tracing in the vacuum Kerr 
space-time could be used for calculating the transfer functions. When computing the 
transfer functions, it was supposed that the matter in the disc rotates on stable circular 
(free) orbits above the marginally stable orbit. The matter below this orbit is freely falling 
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and has the same energy and angular momentum as the matter which is on the marginally 
stable orbit. 

The observer is placed in the direction (/? = 7r/2. The black hole rotates counter- 
clockwise. All six functions are present in these tables. 

Tables are calculated for these values of the black-hole horizon: 

- KBHtablesOO.fits: 1.00, 1.05, 1.10, 1.15, 1.90, 1.95, 2.00 (21 elements), 

- KBHtables50.fits: 1.00, 1.10, 1.20, 1.90, 2.00 (11 elements), 

- KBHtables99.f its: 1.05 (1 element), 

and for these values of the observer's inclination: 

- KBHtablesOO.fits: 0.1, 1, 5, 10, 15, 20, . . . , 80, 85, 89 (20 elements), 

- KBHtables50.fits: 0.1, 1, 10, 20, 80, 89 (11 elements), 

- KBHtables99.f its: 0.1, 1, 5, 10, 15, 20, ... , 80, 85, 89 (20 elements). 

The radii and azimuths at which the functions are evaluated are same for all three 
tables: 

- radii r — are exponentially increasing from to 999 (150 elements), 

- values of the azimuthal angle (/9k are equidistantly spread from to 27r radians with 
a much denser cover "behind" the black hole, i.e. near (/jk = 1.57r (because some of the 
functions are changing heavily in this area for higher inclination angles, > 70°) (200 
elements) . 

B.2 Tables in KBHlineNN . f its 

Pre-calculated functions dF{g) = dg F{g) defined in the eq. (|2.16j) are stored in FITS files 
KBHlineNN . fits. These functions are used by all axisymmetric models. They are stored 
as binary extensions and they are parametrized by the value of the observer inclination 
angle Oq and the horizon of the black hole rt. Each extension provides values for different 
radii, which are given in terms of r — r^, and for different g(-factors. Values of the (^-factor, 
radius r — rh, horizon r^, and inclination 9o, at which the functions are evaluated, are 
defined as vectors at the beginning of the FITS file. 
The definition of the file KBHlineNN . fits: 

0. All of the extensions defined below are binary. 

1. The first extension contains one row with three columns that define bins in the 
5'- factor: 

- integer in the first column defines the width of the bins (0 - constant, 1 - expo- 
nentially growing), 

- real number in the second column defines the lower boundary of the first bin 
(minimum of the ^-factor), 

- real number in the third column defines the upper boundary of the last bin (max- 
imum of the (7- factor). 

2. The second extension contains a vector of the values of the radius relative to the 
horizon r — rh in GM/(P. 

3. The third extension contains a vector of the horizon values in GM/(? (1.00 < rh < 
2.00). 

4. The fourth extension contains a vector of the values of the observer's inclination 
angle Oq in degrees (0° < 60 < 90°, 0° ~ axis, 90° - equatorial plane). 

5. All the previous vectors have to have values sorted in an increasing order. 
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6. In the following extensions the functions dF{g) are defined, each extension is for 
a particular value of and Oq. The values of and 60 are changing with each 
extension in the same order as in tables in the KBHtablesNN.f its file (see the 
previous section, point 7.). Each extension has one column. 

7. Each row corresponds to a particular value of r — rh (see 2. above). 

8. Each element corresponding to a particular column and row is a vector. Each element 
of this vector corresponds to a value of the function for a particular bin in the g- 
factor. This bin can be calculated from number of elements of the vector and data 
from the first extension (see 1. above). 

We have pre-calculated several sets of tables for different limb darkening/brightening 
laws and with different resolutions. All of them were calculated from tables in the 
KBHtablesOO . f its file (see the previous section for details) and therefore these tables 
are calculated for the same values of the black-hole horizon and observer's inclination. All 
of these tables have equidistant bins in the g-factor which fall in the interval (0.001, 1.7). 
Several sets of tables are available: 

- KBHlineOO . f its for isotropic emission, see eq. (|!-{.5|) . 

- KBHlineOl . f its for Laor's limb darkening, see eq. H3.6() . 

- KBHline02 . f its for Haardt's limb brightening, see eq. (|H.7|) . 

All of these tables have 300 bins in the g-factor and 500 values of the radius r — which 
are exponentially increasing from to 999. We have produced also tables with a lower 
resolution ~ KBHlineSO . f its, KBHlineSl . f its, and KBHline52 . f its with 200 bins in 
the g'-factor and 300 values of the radius. 

B.3 Lamp-post tables in lamp, fits 

This file contains pre-calculated values of the functions needed for the lamp-post model. 
It is supposed that a primary source of emission is placed on the axis at a height h above 
the Kerr black hole. The matter in the disc rotates on stable circular (free) orbits above 
the marginally stable orbit and it is freely falling below this orbit where it has the same 
energy and angular momentum as the matter which is on the marginally stable orbit. It 
is assumed that the corona between the source and the disc is optically thin, therefore 
ray-tracing in the vacuum Kerr space-time could be used for computing the functions. 

There are five functions stored in the lamp. fits file as binary extensions. They are 
parametrized by the value of the horizon of the black hole r^, and height h, which are 
defined as vectors at the beginning of the FITS file. Currently only tables for = 1-05 (i.e. 
a = 0.9987492) and h = 2, 3, 4, 5, 6, 8, 10, 12, 15, 20, 30, 50, 75 and 100 are available. 
The functions included are: 

- angle of emission in degrees - the angle under which a photon is emitted from a primary 

source placed at a height h on the axis above the black hole measured by a local 
stationary observer (0° - a photon is emitted downwards, 180° - a photon is emitted 
upwards), 

- radius - the radius in GM/c? at which a photon strikes the disc, 

- ^-factor - the ratio of the energy of a photon hitting the disc to the energy of the same 

photon when emitted from a primary source, 

- cosine of the incident angle - an absolute value of the cosine of the local incident angle 

between the incident light ray and local disc normal, 
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- azimuthal incident angle in radians - the angle between the projection of the three- mo- 
mentum of the incident photon into the disc (in the local rest frame co-moving with 
the disc) and the radial tetrad vector. 

For mathematical formulae defining the functions see eqs. (|2.29|) ~ (|2.3H1 in Section f2.41 
The definition of the file lamp. fits: 

0. All of the extensions defined below are binary. 

1. The first extension contains a vector of the horizon values in GM/c^, though cur- 
rently only FITS files with tables for one value of the black-hole horizon are accepted 
(1.00 < Th < 2.00). 

2. The second extension contains a vector of the values of heights /i of a primary source 
in GM/c^. 

3. In the following extensions the functions are defined, each extension is for a partic- 
ular value of and h. The values of rh and h are changing with each extension in 
the following order: 

rh[l] X h[ll 
rh[l] X /i[2], 
rh[l] X /i[3], 



rh[2] X h[l], 
rh[2] X h[2l 
rh[2] X /i[3], 



4. Each of these extensions has five columns. In each column, a particular function is 
stored - the angle of emission, radius, (^-factor, cosine of the local incident angle and 
azimuthal incident angle, respectively. The extensions may have a different number 
of rows. 



B.4 Coefficient of reflection in f luorescent_line . f its 

Values of the coefficient of refiection f{fj,[,fj,e) for a fluorescent line are stored for different 
incident and refie ction angles in this file. For details on the model of scattering used for 
computations see iMatt. Perola fc Pirol (|l99lh . It is assumed that the incident radiation 



is a power law with the photon index F = 1.7. The coefficient does not change its angu- 
lar dependences for ot her photon indices, only its normalization changes (see Fig. 14 in 
George fc Fabiaiil 1991 ). The FITS file consists of three binary extensions: 



the first extension contains absolute values of the cosine of incident angles, 

the second extension contains values of the cosine of reflection angles, 

the third extension contains one column with vector elements, here values of the 
coefficient of reflection are stored for different incident angles (rows) and for different 
refiection angles (elements of a vector). 
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B.5 Tables in ref spectra . fits 

The function f{Ei]Hi,iie) which gives dependence of a locally emitted spectrum on the 
angle of incidence and angle of emission is stored in this FITS file. The emission is induced 
by a power-law incident radiation. Values of this funct ion were computed by the Monte 
Carlo simulations of Compton scattering, for details see lMatt. Perola Pir 3 (|l99lh . The 



reflected radiation depends on the photon index T of the incident radiation. There are 
several binary extensions in this fits file: 

- the first extension contains energy values in keV where the function /{Ey, ^i,/Xe) is 
computed, currently the interval from 2 to 300 keV is covered, 

- the second extension contains the absolute values of the cosine of the incident angles, 

- the third extension contains the values of the cosine of the emission angles, 

- the fourth extension contains the values of the photon indices T of the incident power 
law, currently tables for T = 1.5, 1.6, . . . , 2.9 and 3.0 are computed, 

- in the following extensions the function /(£"!; /ig) is defined, each extension is for 
a particular value of F; here values of the function are stored as a vector for different 
incident angles (rows) and for different angles of emission (columns), each element 
of this vector corresponds to a value of the function for a certain value of energy. 



Appendix C 

Description of the integration 
routines 

Here we describe the technical details about the integration routines, which act as a 
common driver performing the ray-tracing for various models of the local emission. The 
description of non-axisymmetric and axisymmetric versions are both provided. An appro- 
priate choice depends on the form of intrinsic emissivity. Obviously, non-axisymmetric 
tasks are computationally more demanding. 

C.l Non-axisymmetric integration routine ide 

This subroutine integrates the local emission and local Stokes parameters for (partially) 
polarized emission of the accretion disc near a rotating (Kerr) black hole (characterized 
by the angTilar momentum a) for an observer with an inclination angle Oq. The subroutine 
has to be called with ten parameters: 

ide (ear , ne , nt , far , qar , uar , var , ide_par am , emissivity , neJ-oc) 

ear - real array of energy bins (same as ear for local models in XSPEC), 

ne - integer, number of energy bins (same as ne for local models in XSPEc), 

nt - integer, number of grid points in time (nt = 1 means stationary model), 

far (ne ,nt) - real array of photon flux per bin (same as photar for local models in XSPEC 

but with the time resolution), 

qar (ne , nt ) - real array of the Stokes parameter Q divided by the energy, 

uar (ne , nt ) - real array of the Stokes parameter U divided by the energy, 

var(ne,nt) - real array of the Stokes parameter V divided by the energy, 

ide_param - twenty more parameters needed for the integration (explained below), 

emissivity - name of the external emissivity subroutine, where the local emission of the 
disc is defined (explained in detail below), 

ne_loc - number of points (in energies) where local photon flux (per keV) in the emissivity 
subroutine is defined. 

The description of the ide_param parameters follows: 
ide_param(l) - a/M - the black-hole angular momentum (0 < a/M < 1), 



C.l. Non-axisymmetric integration routine ide 
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ide_parain(2) - theta_o - the observer inclination in degrees (0° - pole, 90° - equatorial 
plane) , 

ide_param(3) - rin-rh - the inner edge of the non-zero disc emissivity relative to the 
black- hole horizon (in GM/c^), 

ide_paraiii(4) ~ ms - determines whether wc also integrate emission below the marginally 
stable orbit; if its value is set to zero and the inner radius of the disc is below the 
marginally stable orbit then the emission below this orbit is taken into account, if 
set to unity it is not, 

ide_paraiii(5) - rout-rh - the outer edge of the non-zero disc emissivity relative to the 
black- hole horizon (in GM/c^), 

ide_param(6) - phi - the position angle of the axial sector of the disc with non-zero 
emissivity in degrees, 

ide_parain(7) - dphi - the inner angle of the axial sector of the disc with non-zero 
emissivity in degrees (dphi < 360°), 

ide_param(8) - nrad - the radial resolution of the grid, 

ide_param(9) - division - the switch for the spacing of the radial grid (0 - equidistant, 

1 - exponential), 
ide_param(10) - nphi - the axial resolution of the grid, 

ide_param(ll) - smooth - the switch for performing simple smoothing (0 - no, 1 - yes), 

ide_param(12) - normal — the switch for normalizing of the final spectrum, 
if = - total flux is unity (usually used for the line), 

if > - flux is unity at the energy = normal keV (usually used for the continuum), 
if < - final spectrum is not normalized, 

ide_param(13) - zshif t - the overall redshift of the object, 

ide_param(14) - ntable - tables to be used, it defines a double-digit number NN in the 
name of the FITS file KBHtablesNN.f its containing the tables (0 < ntable < 99), 

ide_param(15) - edivision- the switch for spacing the grid in local energies (0 - equidis- 
tant, 1 - exponential), 

ide_param(16) - periodic - if set to unity then local emissivity is periodic if set to zero 
it is not (need not to be set if nt = 1), 

ide_param(17) - dt - the time step (need not to be set if nt = 1), 

ide_param(18) - polar - whether the change of the polarization angle and/or azimuthal 

emission angle will be read from FITS tables (0 - no, 1 - yes), 
ide_param(19) - rO-rh and 

ide_param(20) - phiO - in dynamical computations the initial time will be set to the time 
when photons emitted from the point [rO, phiO] on the disc (in the Boyer-Lindquist 
coordinates) reach the observer. 

The subroutine ide needs an external emissivity subroutine in which the local emission 
and local Stokes parameters are defined. This subroutine has twelve parameters: 

emissivity (ear J.OC , ne_loc , nt , f ar_loc , qar_loc , uar_loc , var_loc , r , phi , cosine , 
phiphoton , f irst_emis) 

ear_loc(0:ne_loc) - real array of the local energies where the local photon flux far_loc 
is defined, with special meaning of ear_loc(0) - if its value is larger than zero then 
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the local emissivity consists of two energy regions where the flux is non-zero; the 
flux between these regions is zero and ear_loc(0) defines the number of points in 
local energies with the zero local flux, 
ne_loc - integer, the number of points (in energies) where the local photon flux (per keV) 

is defined, 

nt - integer, the number of grid points in time (nt = 1 means stationary model), 
far_loc(0:ne_loc,nt) - real array of the local photon fiux (per keV) - if the local 
emissivity consists of two separate non-zero regions (i.e. ear_loc(0) > 0) then 
f ar_loc (0 , it) is the index of the last point of the first non-zero local energy region, 

qar_loc(ne_loc,nt) - real array of the local Stokes parameter Q divided by the local 

energy, 

uar_loc(ne_loc,nt) - real array of the local Stokes parameter U divided by the local 
energy, 

var_loc(ne_loc,nt) - real array of the local Stokes parameter V divided by the local 
energy, 

r - the radius in GM/c? where the local photon flux f ar_loc at the local energies ear_loc 

is wanted 

phi - the azimuth (the Boyer-Lindquist coordinate (p) where the local photon flux f ar_loc 

at the local energies ear_loc is wanted, 
cosine - the cosine of the local angle between the emitted ray and local disc normal, 
phiphoton - the angle between the emitted ray projected onto the plane of the disc (in 

the local frame of the moving disc) and the radial component of the local tetrad (in 

radians) , 

f irst_emis - boolean, TRUE if we enter the emissivity subroutine from the subroutine 
ide for the first time, FALSE if we have already been in this subroutine (this is 
convenient if we want to calculate some initial values when we are in the emissivity 
subroutine for the first time, e.g. trajectory of the falling spot). 

C.2 Axisymmetric integration routine idre 

This subroutine integrates the local axisymmetric emission of an accretion disc near a 
rotating (Kerr) black hole (characterized by the angular momentum a) for an observer 
with an inclination angle O^- The subroutine has to be called with eight parameters: 

idre (ear , ne , photar , idre^jaram , cmodel , ne JLoc , ear_loc , f ar_loc) 

ear - real array of energy bins (same as ear for local models in XSPEC), 

ne - integer, the number of energy bins (same as ne for local models in XSPEc), 

photar - real array of the photon flux per bin (same as photar for local models in XSPEc), 

idre_paraiii - ten more parameters needed for the integration (explained below), 

cmodel - 32-byte string with a base name of a FITS flle with tables for axisymmetric 
emission (e.g. "KBHline" for KBHlineNN . fits), 

ne_loc - the number of points (in energies) where the local photon flux (per keV) is 

defined in the emissivity subroutine, 
ear_loc - array of the local energies where the local photon flux f ar_loc is deflned. 



C.2. Axisymmetric integration routine idre 
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f ar_loc - array of the local photon flux (per keV). 

The description of the idre.param parameters follows: 

idre_param(l) - a/M - the black-hole angular momentum (0 < a/M < 1), 

idre_param(2) - theta_o - the observer inclination in degrees (0° - pole, 90° - equatorial 

plane) , 

idre_parain(3) - rin-rh ~ the inner edge of the non-zero disc emissivity relative to the 
black- hole horizon (in GM/c^), 

idre_param(4) - ms - determines whether we also integrate emission below the marginally 
stable orbit; if its value is set to zero and the inner radius of the disc is below the 
marginally stable orbit then the emission below this orbit is taken into account, if 
set to unity it is not, 

idre_parain(5) - rout-rh - the outer edge of the non-zero disc emissivity relative to the 
black-hole horizon (in GM/c^), 

idre_param(6) - smooth - the switch for performing simple smoothing (0 - no, 1 - yes), 
idre_parani(7) - normal — the switch for normalizing the final spectrum, 
if = - total flux is unity (usually used for the line), 

if > - flux is unity at the energy = normal keV (usually used for the continuum), 
if < - final spectrum is not normalized, 

idre_param(8) - zshif t - the overall redshift of the object, 

idre_param(9) - ntable tables to be used, it defines a double-digit number NN in the 
name of the FITS file (e.g. in KBHlineNN.f its) containing the tables (0 < ntable < 
99), 

idre_param(10) - alpha ~ the radial power-law index. 

The subroutine idre does not need any external emissivity subroutine. 



Appendix D 

Atlas of transfer functions 



In this Appendix we provide a graphical representation of the transfer functions that are 
needed for computations of relativistic spectral profiles. For a mathematical definition of 
these functions and details of the adopted notation see Chapter ^ Although graphical 
representation is not necessary in the process of computations and actual data fitting, 
we find it extremely useful and practical for quick order-of-magnitude estimates. Often, 
expected values of various quantities can be seen almost instantly. For example, the mag- 
nitude and the range of the energy shift, the effect of gravitational lensing, the importance 
of the relative time delay etc. can be estimated from these graphs, assuming only basic 
parameters of the light emitting region. 

The following contour graphs of the functions represent a top view of the equatorial 
plane for three different observer's inclination angles (^o = 0.1°, 45° and 85°) in four dif- 
ferent spatial scales (three columns in Figs. ID . itiD .211 and the coarsest scale in Figs. ID.221 - 
ID.23|) . In the finest scale (the first column) we define the radial coordinate in a different 
way to that in the other columns. Specifically, we use r'^ = x'^ -|- y'^ = (r — rh)^. The 
use of r' brings the horizon rh to the origin so that the region just outside the black hole 
is well resolved in these plots. We show figures for two values of an angular momentum 
o of the black hole (rows in Figs. ID.iljD.21|) . particularly for the Schwarzschild black hole 
(a = 0, the horizon is at = 2 and the marginally stable orbit is at r-^as = 6) and for an 
almost extreme Kerr black hole (a = 0.9987, the horizon is at = 1.05 and the marginally 
stable orbit is at rms = 1.198). The clock-wise distortion of the contours visible in the 
second row is due to the frame-dragging near a rapidly rotating Kerr black hole, which is 
clearly visible in the Boyer-Lindquist coordinates. This effect is largely eliminated by the 
transformation from the Boyer-Lindquist to the Kerr ingoing coordinates (the last row of 
panels in Figs. lDlHa2T]) . 

An observer is located to the top of the graphs. The black hole rotates counter- 
clockwise. The values of the functions are encoded by a colour scale, as indicated above 
each panel. The marginally stable orbit is also shown (drawn as a circle) where relevant. 

The atlas shown here is only a part of a larger set that contains figures for more 
values of the observer's inclination and angular momentum of the black hole (including an 
extreme Kerr black hole). An interested reader can ask for the whole atlas by contacting 
the author at |dovcia k@mbo x . troja.mf f . cuni . cz[ 
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Figure D.l: ^(-factor for the inclination 9o = 0.1°. 
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Figure D.2: (^-factor for the inclination Oq = 45°. 
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Figure D.3: (^-factor for the inclination Oq = 85°. 
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Figure D.4: Emission angle for the inclination 9o = 0.1°. 
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Figure D.5: Emission angle for the inclination 9o = 45°. 
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Figure D.6: Emission angle for the inclination 9o = 85°. 
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Figure D.7: Lensing for the inclination 9o = 0.1°. 
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Figure D.8: Lensing for the inclination 9o = 45°. 
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Figure D.9: Lensing for the inclination 9o = 85°. 
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Figure D.IO: Relative time delay for the inclination Oq = 0.1°. 



99 



fi=0: y„^^ Bovflf-Li ndgu i£t coordina te 




^59S7.. J „=.! j05, j-^=l ,1 gs, Bo>B r-Lindi]u isl coord i na Laa 




(i=0,S!9S7, ti,=l JQ5, j-;„^l,JSS, K^rvingoingooofdinaLss 
7J j(i9^2^d^_a&^ 4J^ -Sl4 3H iJ.T JS^ 407 -1 5.7 - 1J.7 -]j5 3H4 




Figure D.ll: Relative time delay for the inclination 0^ = 45°. 
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Figure D.12: Relative time delay for the inclination Oq = 85°. 



101 



fl=0: ■^„m=fi. Bover-Li ndgu ist coordina te 



-laj -90 90 ISO -190 -so 90 133 -ISO -90 90 ISO 




^59S7.. J „=1 j05, j-^=l ,1 gs, BovB r-Lindi]u isl coord i na Laa 




-£i* -at -00 ai 01 HS -.4 -i o i .4 ^ -io -id o lo m 



rtt:0,99S7, j-i,=l i05.. 1 , 1 SS, rr i ngoi ng ooof d i na Lss 



-ISO -90 90 133 -ISO -90 □ 90 ISO -ISO -90 90 IBU 




Figure D.13: Change of the polarization angle for the inclination 9o = 0.1°. 
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Figure D.14: Change of the polarization angle for the inclination 9o = 45°. 
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Figure D.15: Change of the polarization angle for the inclination 9o = 85°. 
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Figure D.16: Azimuthal emission angle for the inclination = 0.1°. 
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Figure D.17: Azimuthal emission angle for the inclination 6o = 45°. 
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Figure D.18: Azimuthal emission angle for the inclination 9o = 85°. 
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Figure D.19: Overall transfer function for the inclination 9o = 0.1°. 
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Figure D.20: Overall transfer function for the inclination Oq = 45°. 
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Figure D.21: Overall transfer function for the inclination Oq = 85°. 
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Figure D.22: 17- factor, emission angle and lensing for the inclinations 9o = 0.1°, 45° and 
85° (from left to right). The strange features in the graph for lensing with the observer 
inclination 0.1° are due to numerical errors which arise when integrating geodesies pass- 
ing close to the axis in spherical coordinates. Notice, though, that the errors are not 
significantly large. 
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Figure D.23: Relative time delay, change of the polarization angle and azimuthal emission 
angle for the inclinations 9o = 0.1°, 45° and 85° (from left to right). 
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